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Abstract 

After giving an outline of the quantization scheme based on the microscopic 
Hopfield model of a dielectric bulk material, we show how the classical phe- 
nomenological Maxwell equations of the electromagnetic field in the presence of 
dielectric matter of given space- and frequency-dependent complex permittivity 
can be transferred to quantum theory. Including in the theory the interaction of 
the medium-assisted field with atomic systems, we present both the minimal- 
coupling Hamiltonian and the multipolar-coupling Hamiltonian in the Coulomb 
gauge. To illustrate the concept, we discuss the input-output relations of ra- 
diation and the transformation of radiation-field quantum states at absorbing 
four-port devices, and the spontaneous decay of an excited atom near the sur- 
face of an absorbing body and in a spherical micro-cavity with intrinsic material 
losses. Finally, we give an extension of the quantization scheme to other media 
such as amplifying media, magnetic media, and nonlinear media. 



*In Coherence and Statistics of Photons and Atoms, edited by J. Pefina (Wiley, New York, 2001), 
p. 1. 
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1 Introduction 

As is already known from classical optics, the use of instruments in optical experi- 
ments needs careful examination with regard to their action on the light under study. 
Fibres, beam splitters, cavities, spectral filters etc. are familiar examples of optical 
instruments, which are typically built up by dielectric bodies. In quantum optics an 
important aspect is the influence of such material bodies on the quantum features 
of light including its interaction with (microscopic) atomic systems. For example, let 
us consider a 50%/50% beam splitter oriented at 45° to an incident light beam. In 
classical optics the beam splitter simply divides the incoming beam into two (apart 
from a phase shift) equal outgoing parts propagating perpendicular to each other, and 
with the same scaling factor the classical noise of the incident field is transferred to 
the two fields in the output channels of the beam splitter. It is intuitively clear that in 
quantum optics the noise of the vacuum in the unused input port of the beam splitter 
introduces additional noise in the two output beams and thus the quantum statistics 
of the output field may differ significantly from that of the input field provided that 
the input field is prepared in a nonclassical state. 

Moreover, when two light beams that are prepared in non-classical states are 
superimposed by the beam splitter, then the outgoing field is prepared in a nonclas- 
sically correlated bipartite state, also called an entangled state. Entanglement as a 
typical quantum coherence feature plays an important role in quantum communica- 
tion. Let us assume that two light beams have been prepared in an entangled state. In 
order to use them, e.g., in quantum teleportation, the beams should be transmitted 
through optical channels such as fibres. Here the crucial point is to what extent the 
entanglement can be preserved during the propagation of the beams, because in any 
real fibre the absorption gives necessarily rise to an entanglement degradation. 

From a more theoretical point of view, a very interesting question is that of the 
Casimir force between material bodies. A crude physical explanation of the Casimir 
force is that the vacuum energies of two regions spatially separated by the bodies 
differ due to the presence of the matter, which in a rough approximation can be 
described simply by a boundary condition on the quantization volume. In fact, this 
simple model does not take account of the dispersive and absorptive properties of the 
matter and fails in the high-frequency limit where the matter becomes transparent. 

These few examples show that it is necessary to include in the theory the presence 
of material bodies when considering the quantized radiation field and its interaction 
with atoms. In principle, material bodies could be included as a part of the matter 
to which the radiation field is coupled and treated microscopically. However, there 
is a class of material bodies whose action can be included in the quantum theory 
exactly, namely dielectric bodies that respond linearly to the electromagnetic field, 
the response being described in terms of a phenomenologically introduced space- 
dependent dielectric permittivity. Such a concept has - similar to classical optics - 
the benefit of being universally valid, because it uses only general physical properties, 
without the need of involved ab initio calculations. 

The quantum theory of radiation in the presence of dielectric matter has been 
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studied over a long period. Quantization of the electromagnetic field in dielectric 
media with assumed real and frequency-independent permittivity has been treated 
extensively [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. In the same context, 
dispersive dielectrics have been considered [16, 17, 18, 19, 20, 21, 22, 23] and at- 
tempts have been made to extend the concepts to nonlinear media [24, 25, 26, 27]. 
However, it is well known that the permittivity is a complex function of frequency 
which has to satisfy the Kramers-Kronig relations, which state that the real part of 
the permittivity (responsible for dispersion) and the imaginary part (responsible for 
absorption) are necessarily connected with each other. A consequence of the existence 
of the imaginary part of the permittivity is that the commonly used mode expansion 
of the (macroscopic) electromagnetic field fails, at least in frequency intervals where 
the absorption cannot be disregarded. Obviously, an expansion of the field in terms of 
damped (non-orthogonal) waves would not be complete. From statistical mechanics 
it is clear that dissipation is unavoidably connected with the appearance of a random 
force which gives rise to an additional noise source of the electromagnetic field. Hence, 
any quantum theory that is based on the assumption of a real permittivity can only 
be valid for narrow-bandwidth fields far from medium resonances where absorption 
can safely be disregarded. 

For the last years there has been an increasing number of articles deal- 
ing with the problem of the formulation of quantum electrodynamics in dielec- 
tric media of given complex permittivity satisfying the Kramers-Kronig relations 
[28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45]. A systematic 
and quantum-theoretically consistent approach to the problem of the quantization 
of the radiation field in absorbing bulk dielectrics is given in [30] on the basis of 
the microscopic Hopfield model of a dielectric [46]. It is based on an explicit Fano- 
type diagonalization [47] of a Hamiltonian consisting of the electromagnetic field, a 
(harmonic-oscillator) polarization field representing the dielectric matter, and a con- 
tinuous set of (harmonic-oscillator) reservoir variables accounting for absorption. The 
resulting expression for the vector potential can be written in terms of the Green ten- 
sor of the classical scattering problem, which, in fact, makes it possible to perform the 
quantization of the electromagnetic field in the presence of arbitrary dielectric bodies 
of (phenomenologically) given permittivities [33, 34, 48], without referring to specific 
microscopic models of the bodies, which were hard to establish for general systems. 
The concept is based on a source-quantity representation of the electromagnetic field, 
in which the electromagnetic-field operators are expressed in terms of a continuous 
set of fundamental bosonic fields via the Green tensor of the classical problem. 

Let us give a brief guide to the topics covered. After giving an outline of the 
quantization scheme based on the microscopic Hopfield model of a dielectric bulk 
material (Sec. 2), we show how the classical phenomenological Maxwell equations 
of the electromagnetic field in the presence of dielectric matter of given space- and 
frequency-dependent permittivity can be transferred to quantum theory (Sec. 3). 
For this purpose we first summarize the basic properties of the classical Maxwell 
equations and express the electromagnetic field in terms of the Green tensor and a 
continuous set of appropriately chosen dynamical field variables (Sec. 3.1). We then 
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perform the quantization by identifying the dynamical field variables with bosonic 
fields associated with the elementary excitations of the composed system (Sec. 3.2). 
Having quantized the electromagnetic field, the question arises of how to include in 
the theory the interaction of the medium-assisted field with atomic systems (Sec. 4). 
In order to answer it, we present both the minimal-coupling Hamiltonian (Sec. 4.1) 
and the multipolar-coupling Hamiltonian (Sec. 4.2) in the Coulomb gauge. To illus- 
trate the basic theoretical concept, we discuss a number of applications such as the 
input-output relations of radiation (Sec. 5.1) and the transformation of radiation- 
field quantum states (Sec. 5.2) at absorbing four-port devices, and the spontaneous 
decay of an excited atom near the (planar) surface of an absorbing body (Sec. 6.2) 
and in a (spherical) micro-cavity with intrinsic material losses (Sees. 6.3 and 6.4). For 
the sake of transparency, we perform the calculations for isotropic media and outline 
the extension to other media in a separate section at the end (Sec. 7). 



2 Hopfield model and Fano diagonalization 

Following the quantization scheme in [30], we consider a Hopfield model [46] of a 
bulk dielectric in which N harmonic oscillator fields describing the polarization of 
the dielectric medium are linearly coupled to a continuum of harmonic oscillators 
standing for a reservoir. Such a model leads to an energy flow essentially only in one 
direction, namely from the medium to the reservoir where it "disappears", hence it 
is absorbed. The overall system that consists of the radiation, the dielectric-medium 
polarization, the reservoir, and couplings between them may be regarded as being a 
Hamiltonian system whose Lagrangian reads 

L — J d 3 r C = J d 3 r (£ rad + £ mat + £ int ), (2.1) 

where 

£ rad = \e a [(A + VU) 2 - c 2 (V x A) 2 ] (2.2) 
(U, scalar potential; A, vector potential; c~ 2 = e fio), and 



N 

•-mat 



X;WX?- W ?X?)+ /dc^[X 2 M-c 2 X 2 M], (2.3) 
i=i i 



N ~ N 

Ant = " ai ( A ^ + f/VX ~ dLJ *M Yl 

1=1 n 1=1 



(2.4) 



Here £ ra d and £ mat are respectively the free Lagrangian densities of the radiation 
field 1 and the matter fields [i.e., the medium oscillator fields X; and the reservoir 
oscillator fields X(c<j) with density //], and C int is the interacting part, where the 
medium-field coupling constants a« play the role of (electric) polarizabilities, and 

1 Note that the vector potential fulfills the Coulomb-gauge condition VA = 0. 
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the medium-reservoir coupling constants Vi{uS) are assumed to be square-integrable 
functions of uj. Introducing the canonical momenta 



^ <9£ • 
n = -j- = e A, 
OA 

^— = /iXj - a { A, 



OX 



QM 



dC 



N 



dX(uj) 



(2.5) 
(2.6) 
(2.7) 



it is not difficult to perform the Legendre transformation and construct the Hamilto- 
nian H = H va ^ + H mSut + H [nt of the overall system. 

Since the dielectric medium is assumed to be infinitely extended, it is convenient 
to go to the reciprocal space 2 , 



A(r)^A(k) = J> A (k)e A (k), 



(2.i 



A=l 



X,(r) -X^k) =X i ||(k)K + J]X iA (k)e A (k) 



(2.9) 



A=l 



[/« = k/|k|, e A (k) _L k], and n(k), Qi(k), X(k, uj), Q(k, uj) accordingly, and to intro- 
duce, with regard to quantization, the new variables 



with 



and 



a A (k) 
Mk) 
b x (k,uj) 



_£o_ 

2hkc 



~ k 2 = k 2 + J2 



2huJi 
2Tiuj 



N 9 



kcA x (k) + -n A (k) 

QiXixQa) + - Qix(k) 
A* 

-iujX x (k, uj) + — Q A (k, a;) 
A* 



~ 2 2 , 

Ui = UJ, + 







/z 2 



ft,, (k, w) 



2ftu)j 

A* 
2/iw 



a) i ||X i ||(k) + -Q i ||(k) 
A 4 



— ia;X|| (k, uj) H — (k, a;) 
A* 



(2.10) 
(2.11) 
(2.12) 

(2.13) 

(2.14) 
(2.15) 



2 The spatial Fourier transform F(k) of a function F(r) is defined according to the relation 
F(t) = (2tt)- 3 / 2 / d 3 k F(k)e 4kr . 
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[£|=of + af/( ) ueo)]. 

The fields are now quantized in the familiar way by conversion of the complex 
amplitudes a A (k), 6j A (||)(k), and &A(||)(k, u) into bosonic annihilation operators a A (k), 
6jA(||)(k), and bu\\\(k,uj), and the Hamiltonian can be expressed in terms of the anni- 
hilation and creation operators as follows: 



H = H 1 - + H majk , = H rad + H^ at + H± t , 

{oo N 
I dco hco b\ (k, W )S|| (k, to) + fit,, (k)6i|| (k) 

1=1 

7 N r 

+ dcoJ2 (^||(k)+&i||(-k)) (b\(-k,u>)+b\\Qt,uj) 



(2.16) 



H- 



i E [d 3 kh~kca[(k)a x (k), 

A=l ^ 



JV 



(2.17) 



(2.18) 



4it = E / d3k | J duJ huJ ^( k > ^ k > ") + E ^ ^(k)fca(k) 

A=l L Q 1=1 

+ /d^E I^H (&L(k)+^(-k)) (Sl(-k,a;)+S A (k,a;)) 

n i=l L 



+EI^H^H (&! A (k)+S iA (-k)j (&] A (-k)+6, A (k; 



(2.19) 



2 N 



H: 



_L 
int 



EE / d3k 5^(*) [4(-k) + a A (k)] [ftt^k) + Sa(-k)] , (2.20) 

A=l i=l 

where Vi{u) = (vi(u) / fai)' 1 / 2 and Ai(k) = [(iJia 2 ) / (fice k)} 1/2 . The bilinear Hamil- 
tonian can be diagonalized (separately for the transverse and longitudinal parts) by 
applying a Fano-type technique [47] such that 3 



H, 



mat 



J d 3 k J duo hco s[(k,w)S||(k,w), 
o 



2 

# x = E / d3 k f dcohcoCl(k,co)C x (k,co). 
a=i ^ * 



(2.21) 



(2.22) 



3 A different technique for diagonalization is utilized in [41, 42], where path- integral quantization 
of the Lagrangian (2.1) is performed. 
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Since the derivation of the formulas for expressing the old bosonic operators frj||(k), 
( k )'^ll ( k > w )> ( k > w ) and 6 *(k, 4 ( k , M k ), 6- A (k), S A (k, w), 6 A (k, u) in terms 
of the new bosonic operators B\\(k.,u), B^(k,u) and C\(k,u), C\(k,cu), respectively, 
is rather lengthy, we renounce it here and refer the reader to [30, 49]. 

The vector potential and the transverse part of the medium polarization can 
then be expressed in terms of the (polariton-like) operators C\(k,u), C\(k,u) as 
[30, 34, 48] 

2 / 00 

A(k) = - Ve A (k)W— f du uy/eAu) G(k, u)C x (K u) + H.c. (2.23) 

V TT£ J 

A-l q 

and 

2 N 



2 2 fh — 



A=l i=l A=l 



X 



oo 

y dcj^ 2 [e(cj)-l] v^eiH G(k, u)C x (k, u) +H.c. 



P N (k), (2.24) 



2 / 

P N (k) = ^e A (k)^/^ Je^iy/^(ujC x (k,u) + H.c.. (2.25) 

A=l 

Here, s{uj) —e^(cu) +iei(u>) is the complex (model) permittivity of the medium, and 



c 2 



GM = - J>e( U )-W (2 ' 26) 

is the Green function of the classical Maxwell equations with that permittivity. From 
Eq. (2.24) [together with Eq. (2.25)] it is seen that the (transverse) polarization of 
the medium consists of two qualitatively different terms. Obviously, the first term 
is the induced polarization, whose frequency components are given by the frequency 
components of the electric field strength of the radiation multiplied by eo[e(u) — 1], 
and the second term is the (noise) polarization, i.e., the fluctuating component of the 
polarization that is associated with absorption. 

Recalling that the longitudinal part of the electric field operator is given by 

eoEii(k) = — Q!jXj||(k)K and recalling the transversality of the displacement op- 
erator D(k), a relation between Ey and a longitudinal (noise) polarization defined 
analogously to Eq. (2.25) [by replacing C\(k,u) with £||(k, cu)] can be derived. 

In summary, the quantized electromagnetic field can be expressed, via a source- 
quantity representation with the classical Green function, in terms of the permittivity 
and a continuum of harmonic oscillators. It is worth noting that in this formulation 
there is no explicit hint at the underlying microscopic model. Thus, it seems quite 
natural to generalize the theory to arbitrary dielectric matter of given space- and 
frequency-dependent permittivity by transferring the classical source-quantity repre- 
sentation of the electromagnetic field directly to quantum theory. 



3 THE MEDIUM-ASSISTED MAXWELL FIELD 



8 



3 The medium-assisted Maxwell field 

From now we will not refer to one or the other microscopic model of the dielectric 
media. Instead we will start from the familiar phenomenological Maxwell equations, 
assuming that the permittivity is known, e.g., from measurements. In order to allow 
for an arbitrary formation of different (non-moving) dielectric bodies in space, we 
will assume that the permittivity varies with space. For the sake of transparency we 
will disregard the tensor character of the permittivity, restricting our attention to 
isotropic media (for the extension to anisotropic media, see Sec. 7.2). 



3.1 Classical basic equations 

Let us first briefly outline the classical theory and bring it in a form suitable for 
quantization. The phenomenological Maxwell equations of the electromagnetic field 
in the presence of dielectric bodies but without additional charge and current densities 
read 

VB(r) = 0, (3.1) 

V x E(r) + B(r) = 0, (3.2) 

VD(r) = 0, (3.3) 

V x H(r) - D(r) = 0, (3.4) 

where the displacement field D is related to the electric field E and the polarization 
field P according to 

D(r)= £o E(r) + P(r), (3.5) 
and for nonmagnetic matter it may be assumed that 

H(r) = -B(r) (3.6) 

(for the extension to magnetic matter, see Sec. 7.3). 

Let us consider arbitrarily inhomogeneous (isotropic) media and assume that the 
polarization responds linearly and locally to the electric field. In this case, the most 
general relation between the polarization and the electric field which is in agreement 
with the causality principle and the dissipation-fluctuation theorem is 

POO 

P(r,t)=e / dr X (r,r)E(r,t-r) + P N (r,t), (3.7) 
Jo 

where x( r ) T ) is the dielectric susceptibility as a function of space and time, and Pn 
is the (noise) polarization associated with absorption. 

Substitution of this expression into Eq. (3.5) together with Fourier transformation 
converts this equation to 



D(r, uj) = e e(r, w)E(r, u) + P N (r, u), 



(3.8) 
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and thus 

P(r, u) = e [e(r, u) - 1] E(r, u) + P N (r, w), (3.9) 

where 

poo 

e(r,u) = l+ / dr X (r,r)e^ (3.10) 
./o 

is the (relative) permittivity, and the Maxwell equations (3.1) - (3.4) read in the 
Fourier domain as 4 

VB(r,cj) = 0, (3.11) 
VxE(r,w)-iwB(r,w)=0, (3.12) 
£ Ve(r,cu)E(r,cu) = p N (r,cu), (3.13) 

V x B(r, u) + i^e(r, w)E(r, w) = // j N (r, w). (3.14) 

Here we have introduced the (noise) charge density 

P N (r,w) = -VP N (r lW ) (3.15) 

and the (noise) current density 

j N (r, u) = -^P N (r, uj), (3.16) 

which obey the continuity equation 

Vj N (r,u;)-^p N (r,u;) = 0. (3.17) 

According to Eq. (3.10), the permittivity e(r,iv) is a complex function of fre- 
quency, 

e(r,u) =e R (r,u)+iei(r,u). (3.18) 

The real and imaginary parts, which are responsible for dispersion and absorption, 
respectively, are uniquely related to each other through the Kramers-Kronig relations 

£R (r, w )-l = ^/da/^, (3.19) 
e^M = -- fdu,' 6 *^- 1 (3.20) 

7T J Uj' — UJ 

(V, principal value). Further, e(r,u>) as a function of complex uj satisfies the relation 

e(r,-u*) =e*(v,cu) (3.21) 

and is holomorphic in the upper complex half-plane without zeros. In particular, it 
approaches unity in the high-frequency limit, i.e., e(r, uj) — > 1 if \uj\ — > oo [50, 51]. 

4 Here and in the following the Fourier transform F[u>) of a real function F(i) is defined according 
to the relation F(t) = (t) + (t) , where (i) = J °°dw F{oj)e- iut and F^ (t) = [*"(+) (*)]*. 
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The Maxwell equations (3.12) and (3.14) imply that E(r, co) obeys the partial 
differential equation 

V x V x E(r,w) - — e(r,cu)E(r,cj) = ^/i j N (r,^), (3.22) 
the solution of which can be represented in the form 

E(r, co) = zufio J dV G(r, r', co)j_ N (r>, co), (3.23) 
where the Green tensor G(r, r', a;) has to be determined from the equation 

CO 2 

V x V x G(r, r', co) - — e(r, w)G(r, r', w) = <5(r - r') (3.24) 



together with the boundary condition at infinity. In Cartesian coordinates, Eq.(3.24) 
reads 



co 2 



(dld r k - 5 lk A r ) - S ik —e(r,u;)\G kj (r,r',uj) = 5 l3 5{v - r') (3.25) 

{d[ —d/dxi), where over repeated vector-component indices is summed. The Green 
tensor has the properties that 

G; j (ry,u) = G ij (ry,-u*), (3.26) 

G'„(r',r,cc;) = G'^(r,r / ,^), (3.27) 

and 

J d 3 s ^ £l (s, co) G lk (r, s, u)G* jk (r', s, uo) = Im G tJ (r, r', cu). (3.28) 

The property (3.26) is a direct consequence of the corresponding relation (3.21) for 
the permittivity. The reciprocity relation (3.27) and the integral relation (3.28) are 
proven in Appendix A. 

The Fourier components of the magnetic induction, B(r, to), and the displacement 
field, D(r, co), are directly related to the Fourier components of the electric field, 
E(r,co), 

B(r, co) = (ico)- 1 ^ x E(r, co), (3.29) 

D(r, co) = (iiquo 2 )- 1 ^ x V x E(r, co) (3.30) 

[see Eqs. (3.12), (3.8), (3.16), and (3.22)], and E(r, co) is determined, according to 
Eq. (3.23), by j N (r, co). The continuous set of (complex) fields j N (r, co) [or equiva- 
lently, P N (r, co)] can therefore be regarded as playing the role of the set of dynamical 
variables of the overall system composed of the electromagnetic field and the medium 
(including the dissipative system). For the following it is convenient to split off some 
factor from P N (r, co) and to define the fundamental dynamical variables f(r, co) ac- 
cording to 

[he 

P N (r,u) = id —^-€i(r,uj)f(r,uj). (3.31) 
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3.2 Field quantization 

The transition from classical to quantum theory now consists in the replacement of 
the classical fields f(r,u;) and f*(r,cj) by the operator-valued bosonic fields f(r,u;) 
and ft(r, a;), respectively which are associated with the elementary excitations of the 
composed system within the framework of linear light-matter interaction. Thus the 
commutation relations are 

[/ fc (r,w),/t (rV)] = 6 kk ,5(r-r')6(uj-uj'), (3.32) 

[f k (r,u)J k ,(r',u')] = 0, (3.33) 
and the Hamiltonian of the composed system is 

/POO 
d 3 r J dwnw^(r,u)i(r,u). (3.34) 

Replacing E(r, uj) [Eq. (3.23)], B(r, u) [Eq. (3.29)], and D(r, u) [Eqs. (3.8), (3.30)] 
by the quantum-mechanical operators, we find, on recalling Eqs. (3.16) and (3.31), 
that 

E{r,u) = u[-^ < ^ fdVv^MG^r'^jf^o;), (3.35) 

V 7T£ & J 

B(r, uj) = (iw) _1 V x E(r, uj), (3.36) 

and 

D(r, uj) = e e(r, w)E(r, uj) + P N (r, uj) 

= (fi cu 2 y l V x V x E(r, uj), (3.37) 

from which the electromagnetic field operators in the Schrodinger picture are obtained 
by integration over uj: 

poo 

E(r) = / dwE(r,o;) + H.c., (3.38) 
Jo 

POO 

B(r) = / dwB(r,w) + H.c, (3.39) 

POO 

D(r) = D-^r) = / dwD(r,w) + H.c. (3.40) 

In this way, the electromagnetic field is expressed in terms of the classical Green tensor 
G(r, r',uj) satisfying the generalized Helmholtz equation (3.24) and the continuum 
of the fundamental bosonic field variables f(r,u;) [and F(r, uj)]. All the information 

5 Here the longitudinal (FH) and transverse (F^) parts of a vector field F are defined by FlK-^fr) 
= Jd 3 r' ^"^ ± - ) (r — r')F(r'), with 5^"^(r) and 5^(r) being respectively the longitudinal and trans- 
verse tensor-valued (5-functions [Eqs. (B.9) and (B.10)]. Note that for bulk material the transverse 
polarization field = D — eoE^, with E and D being respectively given by Eqs. (3.38) and (3.40) 
together with Eq. (3.35) and (3.37) exactly corresponds to Eq. (2.24) together with Eq. (2.25). 
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about the dielectric matter (such as its formation in space and its dispersive and 
absorptive properties) is contained [via the permittivity e(r,u)] in the Green tensor 
of the classical problem. Eqs. (3.38) - (3.40), together with Eqs. (3.35) - (3.37), can 
be considered as the generalization of the familiar mode decomposition. 

A similar formalism which also starts from a causal relation between the polar- 
ization and the electric field strength is developed in [37, 38]. The auxiliary fields 
that are introduced there in order to construct a unitary time evolution in an en- 
larged Hilbert space can be shown to lead essentially to the field variables f(r, ou) 
and f'(r, uf) considered here. Thus, the representation of the electromagnetic field in 
[37, 38] corresponds to the Green function representation in Eqs. (3.35) - (3.40). 

The quantization scheme meets all the basic requirements of quantum electrody- 
namics. So it can be shown by using very general properties of the permittivity and 
the Green tensor that the electric and magnetic fields satisfy the correct (equal-time) 
commutation relations (see Appendix B) 

[E k (r), E t (r')] = = [B k (r), B^r')} , (3.41) 

[e E k (r), = -ihe klm d r J{r-v'). (3.42) 

Obviously, the electromagnetic field operators in the Heisenberg picture satisfy the 
Maxwell equations (3.1) - (3.4), with the time derivative of any operator Q being 
given by 

Q={ihY l [Q,H], (3.43) 

where H is the Hamiltonian (3.34). 

Let us briefly comment on the (zero-temperature) statistical implications of the 
quantization scheme. The vacuum expectation value of E(r, u) is obviously zero 
whereas the fluctuation of E(r, u) is not. From Eq. (3.35) together with the com- 
mutation relations (3.32) and (3.33) we derive, on using the integral relation (3.28), 

<0|4(r,w)4V,a/)|0) = ^ImG H (r,r»^- W '). (3.44) 

Equation (3.44) reveals that the fluctuation of the electromagnetic field is deter- 
mined by the imaginary part of the Green tensor - a result that is consistent with 
the dissipation-fluctuation theorem 6 [52]. Thus, the quantization scheme respects 
both the basic requirements of quantum theory (in terms of the correct commutation 
relations) and statistical physics (in terms of the dissipation-fluctuation theorem). 

So far we have considered the electromagnetic field strengths. Instead, scalar (ip) 
and vector (A) potentials can be introduced and expressed in terms of the fundamen- 
tal bosonic fields f(r, u) and F(r, u). In particular, the potentials in the Coulomb 
gauge are defined by 

-V0(r) = El'(r), (3.45) 

6 Note that the Green tensor plays the role of the response function of the electromagnetic field 
to an external perturbation. 
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A(r) 



du;A(r, uj) + H.c. , 



where 7 



A(r,w) = M _1 E (r,uj). 
The canonically conjugated momentum field with respect to A(r) is 



POO 

n(r) = — ie / dujuj A(r, uj) + H.c. , 



(3.46) 
(3.47) 

(3.48) 



and it is not difficult to verify that n = — eoE- 1 , V x A = B, and —A — Vy3 = E. In 
addition, A and II satisfy the well-known commutation relations (Appendix B) 



[A k (r),A k ,(r')] = 0= [n fc (r),n fc ,(r')], 
[i fc (r),H fe ,(r')] =<R^(r-rO. 



(3.49) 
(3.50) 



3.2.1 One-dimensional systems 



Let us illustrate the main features of the concept for linearly polarized radiation 
propagating in the redirection, which effectively reduces the system to one spatial 
dimension [A(r) — > A y (x) = A{x), n(r)— > fl y (x) = TL(x), f(r, u) — >f(x,uj)\. According 
to Eqs. (3.46) and (3.47) together with Eq. (3.35), the operator of the vector potential 
is 

A(x) = ^~^jJ^ du J Ax ^ e^x' ,uj) G(x,x' ,uj)jXx } uj) + R.c. , (3.51) 

and IT accordingly (.4., normalization area perpendicular to the a; direction). Here the 
Green function G(x, x', uj) satisfies the equation 



d uj 2 

■——G(x, x' , uj) e(x, uj)G(x, x', uj) = 5(x — x'). 

ox 2 c 2 



(3.52) 



In the simplest case when the spatial variation of the permittivity can be disregarded, 
then the solution of Eq. (3.52) that satisfies the boundary conditions at \x\, \x'\ ^oo 
is 



G(x, x' , uj) = — 



.UJ , , 

2i— n(uj) 
. c 



cxp 



i— n(uj)\x - x 
c 



(3.53) 



with n(uj) = ^y e(uj)=nn(uj)+ini(uj) being the complex refractive index of the medium. 

7 Note that for bulk material Eqs. (3.46) and (3.47) together with Eq. (3.35) exactly correspond 
to Eq. (2.23). 
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Substituting the Green function (3.53) into Eq. (3.51), we may rewrite the x'- 
integral to obtain 



4neQCUjn-R{uj)A n(u) 
x [e mn(uj)u)x/c a + (x,u) +e- in ^ u)X / c a^(x,u)] +H.c.|, (3.54) 

where 

a±(x, to) = i^2 ni {u)uj/c e ^(-W c I dx' e- m{uj)uJX ' /c f(±x', to), (3.55) 

J — oo 

and from Eq. (3.32) it follows that 

[a±(x, uj), ai(x', uj')] = e ~M^-A/c _ (3 56 ) 

Obviously, the space-dependent operators a±(x, uj) describe the amplitudes of the 
damped monochromatic waves propagating to the right (subscript +) and left (sub- 
script — ), and from Eq. (3.55) it follows that a±(x,cu) and a±(x', uj) are related by 
(spatial) quantum Langevin equations [30, 33], 

-^a±(x,uj) = Tni(uo)^a±(x,uj) + F±(x,oj), (3.57) 

with 

F±(x, to) = ±i v / 2n l (Lu)iu/c e^'' f(x, u) (3.58) 

being the operator Langevin noise sources. In particular, when (f(x",oj)) = for \x" 
— x\ < \x — x'\ and \x" — x'\ < \x — x'\, then 

(a±(x, uj)) — {a±(x', uj)) exp[— m(u)u\x — x'\/c], ±x =1= x > (3.59) 

for arbitrary (a±(x',uj)). 

Equation (3.54) is the extension of the familiar mode decomposition to absorb- 
ing media. Let us assume that in a frequency interval Au the absorption is suffi- 
ciently small, so that for a chosen (finite) propagation interval \x — x'\ the condition 
ni(uj)uj\x — x'\/c<^ 1 holds. Then the amplitude operators a±(x,uj) can be regarded 
as being independent of x for that propagation distance and satisfying the ordinary 
Bose commutation relations, as is seen from Eqs. (3.55) and (3.56). In the chosen 
frequency and space intervals, Eq. (3.54) exactly reduces to the familiar expression 
obtained by mode expansion. 
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4 Atom— Field Interaction 

The interaction of the quantized electromagnetic field with atoms placed inside a 
dielectric-matter configuration or near dielectric bodies can strongly be influenced 
by the dielectric medium. A well-known example is the dependence of the rate of 
spontaneous decay of an excited atom on the properties of a dielectric environment 
(Sec. 6). In order to study such and related phenomena, the Hamiltonian (3.34) 
must be supplemented with the Hamiltonian of additional charged particles and their 
interaction energy with the medium-assisted electromagnetic field. 

4.1 The minimal-coupling Hamiltonian 

Applying the minimal-coupling scheme, we may write the total Hamiltonian in the 
form 8 

d 3 r J dco huj f f (r, u)f (r, u) + ^ — |p Q - g a A(r a ) 

+ \ J d 3 rp A (r)y3 A (r) + J d 3 r p A (r)y3 M (r), (4.1) 

where r a is the position operator and p a is the canonical momentum operator of 
the ath (nonrelativistic) particle of charge q a and mass m a . The Hamiltonian (4.1) 
consists of four terms. The first term is the energy of the electromagnetic field and the 
medium (including the dissipative system), as introduced in Eq. (3.34). The second 
term is the kinetic energy of the charged particles, and the third term is their Coulomb 
energy, where the corresponding scalar potential <p A is given by 

with 

p A (r) = ^a(r-f«) (4.3) 

a 

being the charge density of the particles. The last term is the Coulomb energy of 
interaction of the particles with the medium. From Eq. (4.1) it follows that the 
interaction Hamiltonian reads 

#int = - — \p a - |g Q A(r a )j A(r a ) + / d 3 r p A (r)y3 M (r). (4.4) 
a m a l J j 

Note that in Eqs. (4.1) and (4.4), the vector potential A and the scalar potential y3 M , 
respectively, must be thought of as being expressed, on using Eq. (3.46) [together 
with Eqs. (3.47) and (3.35)] and Eq. (3.45) [together with Eqs. (3.38) and (3.35)], in 
terms of the fundamental fields f(r,cj) [and f '(r, a;)]. 

8 Here and in the following the subscripts A and M are introduced in order to distinguish between 
atom- and medium-assisted quantities. 



4 ATOM-FIELD INTERACTION 



16 



In a straightforward but somewhat lengthy calculation (for an example, see Ap- 
pendix C) it can be shown (by means of the commutation relations derived in Ap- 
pendix B) that both the operator-valued Maxwell equations 



where 



VB(r) = 0, 

V x E(r) + B(r) = 0, 
VD(r)=p A (r), 

VxH(r)-D(r)=j A (r), 
j A (r) = \ S ^q a [r a 5(r - r Q ) + 5(r - f Q ) i a 



and the operator-valued Newtonian equations of motion 

• 1 



a r a = q a \ E(f Q ) + \ 



p Q - q a A(r a ) 
r n x B(r J - BffJ x r r 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 



(4.10) 
(4.11) 



are fulfilled. In Eqs. (4.6) - (4.8), the (longitudinal part of the) electric field and 
the displacement field now contain [compared with Eqs. (3.38) and (3.40)] additional 
longitudinal components that result from the charge distribution p A (r), i-e., 



E(r) = E M (r) - V^ A (r) 



D(r) = D M (r)-£oV^ A (r) 



dcuE(r,cj) + H.c. 



' POO 

/ dcjD(r,cu) + H.c. 
Jo 



- V^ A (r), 



(4.12) 



£ Vy9 A (r), (4.13) 



with E(r, u) and D(r, uo) being defined by Eqs. (3.35) and (3.37). The Maxwell equa- 
tions (4.5) and (4.7), respectively, simply follow from the definition of B(r) [Eq. (3.39) 
together with Eq. (3.36)] and D(r) [Eq. (4.13) together with Eqs. (3.37) and (4.2)]. 
The Maxwell equations (4.6) and (4.8) are respectively the Heisenberg equations of 
motion of B(r) and D(r) according to Eq. (3.43), with the Hamiltonian being given by 
Eq. (4.1), and the Newtonian equations of motion (4.10) and (4.11) are respectively 
obtained from the Heisenberg equations of motion of r a and p a . 



4.2 The multipolar-coupling Hamiltonian 

In the interaction Hamiltonian (4.4) used in the minimal-coupling scheme the electro- 
magnetic field is expressed in terms of the potentials. With regard to the interaction of 
the electromagnetic field with (localized) atomic systems (atoms, molecules etc.) the 
interaction energy is commonly desired to be treated in terms of the field strengths 
and the atomic polarization and magnetization. This can be achieved by means of a 
unitary transformation. 
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Let us consider an atomic system localized at position r A and introduce the atomic 
polarization 



V q a (r Q - r A ) / dA5[r-r A -A(r a -r A )] , 



so that the charge density (4.3) can be rewritten as 

Pa(t) = q A S(r - r A ) - VP A (r), 

with 



(4.14) 



(4.15) 



(4.16) 



being the total charge of the atomic system. In order to perform the transition from 
the minimal-coupling scheme to the multipolar-coupling scheme, we apply to the 
variables the unitary operator 



U = exp 



- /d 3 rP A (r)A(r) 



(4.17) 



which is known as the Power-Zienau transformation [53, 54, 55, 56]. It is not difficult 
to prove that the following transformation rules are valid: 



f' = Ur = r 



(4.18) 



i>' a = ui> a u ] 

= p a - q a A(r a ) - q a dX X (r a -r A ) x B [r A + A (f Q -r A )] , (4.19) 

Jo 



f'(r,u) = Uf(r,w)(P 



= f(r,u)-U—e l (r,u)^ /dV Pi(r')G*(r', r, w ). (4.20) 

Employing equations (4.18) - (4.20) and using Eqs. (3.28), (3.35), and (B.13), we can 
express the Hamiltonian H in Eq. (4.1) in terms of the new variables r^, = r a , p^, and 
f'(r,u) in order to obtain the multipolar Hamiltonian 

H = Jd 3 rJ^ dujkuji n (r,uj){'(r,uj) 

+ J2^~{^ + qa J o d ^^(^-r A ) x B'[r A +A(r a -r A )] | 
+ \ J d 3 rp A (r)y3 A (r) + ±- J d 3 r P*(r)P*(r) 

+ J d 3 rp A (r)y3Ur) - j d 3 r F»i(r)E^(r), (4.21) 
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where the relations B' = B, (f' M = (p M [cf. Eq. (B.19)], and 

E&(r) = Ei(r) + -F»i(r) (4.22) 

are valid. 

In particular when the charged particles form a neutral atomic system (g A = 0), 
then we may write, on integrating by parts and recalling that = _ Pa(m)/ £ 0' 



/ 



d 3 rp A (r)y3 A (r) = - [ d 3 r P A (r)P A (r), (4.23) 
£o J 



and 



J d 3 rp A (r)^ M (r) = ±- J d 3 r p[(r)p|[ I (r) 

= - / d 3 rP A (r)P„(r) - - I d 3 r P^P^r). (4.24) 

Combining Eqs. (4.21) - (4.24) and taking into account that P^ = Pm [cf. Eqs. (B.14) 
and (B.16)], we may rewrite the multipolar Hamiltonian as 

//»oo 
d 3 r J dLuhu{'\r,u){'(r,Lu) 

+ 5Z^~t{p° + ^/ dAA(r a -r A ) x B' [r A + \ (r a -r A )} j 

+ J_ /d 3 rP A (r)P A (r)+ - /*d 3 rP A (r)P' M (r) 

2^o J £o J 

-- fd s rP A (r)iy ± (r), (4.25) 
£o J 

where 

D'^(r) = D£(r) = ^oE^(r) + Pft(r). (4.26) 
From Eq. (4.25) the interaction Hamiltonian is seen to be 

H int = - I d 3 rP A (r)P^(r) - - I d 3 r P A (r)D' ± (r) 
£o J £o J 

~ E ^ jf dA A {[Pc-ta) x P«l B'[^ + A (r a -r A )} + H.c.} 

+ ^^^J o 1 dXX(r a -r A )xB'[r A + X(r a -r A )]^ . (4.27) 

The first term on the right-hand side in Eq. (4.27) is a contact term between the 
medium polarization and the polarization of the atomic system. The second term 
describes the interaction of the polarization of the atomic system with the transverse 
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part of the overall displacement field [cf. Eq. (4.26)], and the last two terms refer to 
magnetic interactions. 

So far we have transformed the dynamical variables but left unchanged the Hamil- 
tonian. Instead the Hamiltonian can be transformed to obtain the new one 

H = U ] HU. (4.28) 

Obviously, the new Hamiltonian expressed in terms of the new variables formally looks 
like the untransformed minimal-coupling Hamiltonian expressed in terms of the old 
variables. Hence expressing in the new Hamiltonian the new variables in terms of the 
old ones, we arrive at a multipolar-coupling Hamiltonian that (for a neutral atomic 
system) looks like the Hamiltonian given in Eq. (4.25), i.e., 

//»oo 
d 3 r J &ufiwi\r,u)i{r,u) 

+ Y1 2^r{^ Q + <la J dAA (^~ r ^) x B[r4 + A(r a -r A )] j 

+ J_ /d 3 rP A (r)P A (r)+ - fd 3 rP A (r)P M (r) 

2e J e J 

-— /d 3 rPA(r)D- L (r). (4.29) 
£o J 

In fact, the Hamiltonians in Eqs. (4.25) and (4.29) have different meanings. Since the 
expectation value of an observable associated with an operator O must not change, 
the use of 7i necessarily requires a transformation of both the operator [O — > WOU] 
and the state [q—>WqU, with g being the density operator], so that 

ih^(0) =Tf{q[0,H]} =Tr{tfQU[&OU,H\}. (4.30) 

Moreover, in Eq. (4.25) the atomic polarization field couples to the transverse compo- 
nent of the overall displacement field, which also contains the transverse component 
of the atomic polarization, rather then the ordinary transverse displacement field in 
Eq. (4.29). 



5 Input— output coupling 

The quantization scheme developed in Sec. 3 is best suited to study the input-output 
behaviour of optical fields at dielectric devices, because both dispersion and absorp- 
tion are exactly included in the resulting input-output relations, and characteristic 
quantities such as the transmission and reflection coefficients of the setup are ex- 
pressed in terms of its complex refractive- index profile [57]. Input-output relations 
are a very efficient description of the action of macroscopic bodies on radiation. In 
particular, they can advantageously be used to obtain the quantum statistics of the 



5 INPUT-OUTPUT COUPLING 



20 



outgoing light from that of the incoming light, either in terms of radiation- field cor- 
relation functions or directly in terms of the density matrix. In what follows we will 
restrict our attention to four-port devices. The extension of the method to (higher- 
order) multiport-devices is straightforward. 

5.1 Operator input— output relations 

Let us study the problem of propagation (in the x-direction) of quantized radiation 
through a dielectric plate that is the middle part of a three-layered planar structure as 
sketched in Fig. 1. The setup can be characterized by a piecewise constant permittivity 



o o o 
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° efaf 
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• • 
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1 1— 

-1/2 1/2 



Figure 1: Scheme of a three-layered structure: a dielectric plate of thickness I and 
permittivity 82(00) is surrounded by dielectric matter of permittivities £1(0;) (on the 
left-hand side) and £3(0;) (on the right-hand side). 

e(x,u>) = JXs) e », \j(x) = { I f th ^/ < Xj , (5.1) 
3=1 < 

where Sj(uj) is the complex permittivity of the jth domain, and x — > —00, X\ = —1/2, 
X2 — I/2, and £3— > + oo. The Green function G(x, x', uj) satisfies the partial differential 
equation (3.52) from which it follows that it can be decomposed into two parts 9 

3 

G(x,x',uj) =J2 X j( x ) X j( x ') G T k ( x i x 'iu)+R(x,x',L;), (5.2) 
3=1 

where G^ nlk (x, x 1 \ uj) is the Green function of the bulk material as given by Eq. (3.53) 
(with Ej in place of e), and the reflection term R(x,x',uj) is a solution to the homo- 
geneous wave equation, 

3 

R(x, x', uj) = J2 X i( x ) [ C J+( X ', u)e m ^ x/c + Cj4x', w ) e - <n ^ w W c ] , (5.3) 

j=i 

9 For the calculation of Green tensors and a large variety of examples, see [58]. 
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where the coefficients Cj±(x',ou) are to be determined in such a way that continuity 
and differentiability at the surfaces of discontinuity at X — —1/2 and x — l/2 are en- 
sured. Combining Eqs. (3.51) and (5.2) [together with Eq. (5.3)], the vector potential 
A(x) for the jth domain may be represented as, similar to Eq. (3.54), 



A(x) = [ 
Jo 



dui 



h n jR (uj) 



x [e in ' R(w)wx ^^ (5.4) 

[xj-i <x <Xj], where the dependence on x of the amplitude operators a,j±(x,u>) is 
governed by quantum Langevin equations of the type given in Eq. (3.57) together 
with Eq. (3.56). 

The amplitude operators of the incoming fields are given according to Eq. (3.55) 
and satisfy the commutation relations 10 

[a 1+ (x, u), a\ + (x', u')] = e -^{^\x-x'\/c 5( ^ _ ^ ^ 

[a 3 -(x, u), al_(x', u')] = e^M"!*-*'!/^ _ ( 5 . 6 ) 
[a 1+ (x,uj),al_(x',uj')] = (5.7) 

[cf. Eq. (3.56)]. Thus, the incoming fields from the left and right behave like the fields 
in the corresponding bulk dielectrics and may be regarded as independent variables. 
Further, it can be shown that the amplitude operators of the outgoing fields can be 
related to the amplitude operators of the incoming fields and appropriately chosen 
operators g±(w) of the plate as 

where the 2x2-matrices T{uj) and A (a;) are the characteristic transformation and 
absorption matrices of the plate expressed in terms of the thickness and the permit- 
tivity of the plate and the permittivities of the surrounding media, and the operators 
g±{oj) read 



in 2 (u>)u>l/(2c) 



X 



f l/2 

/ dx' [e in2(w)uJX ' /c ±e- in2 ^ )wx ' /c ]f(x f ,to), (5.9) 

J -1/2 



10 Note that the amplitude operators of the field inside the plate (domain 2) and the amplitude 
operators of the outgoing fields are not given according to Eq. (3.55) and thus do not satisfy com- 
mutation relations of this type in general. 
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with 

(for details, see [57]). It is not difficult to prove, on recalling the basic commutation 
relations (3.32) and (3.33), that the commutation relations 

[g ± (u),gi(u')] =8{u-w'), (5.11) 

[g±(u), =0 (5.12) 

are valid. Hence, the operators g±(u>) and g±(w) are respectively annihilation and 
creation operators of bosonic excitations associated with the plate. Obviously, they 
commute with the amplitude operators of the incoming fields a 1+ (x, u) and a 3 _(x, ui) 
and are thus independent variables. It should be mentioned that the output amplitude 
operators a,i-(x,u), x < —1/2, and a 3+ (x, uj), x > 1/2, can easily be obtained from 
a,i-(—l/2,u) and a 3+ (l/2,u), respectively, by means of the corresponding solutions 
of the Langevin equations (3.57). 

Whereas the input amplitude operators a\ + (—l/2,u), cl\ + (—1/2,uj) and 
a^-{l/2,u), a\_{l/2,uj) satisfy ordinary bosonic commutation relations [see Eqs. (5.5) 
- (5.7)], the amplitude operators of the outgoing fields, ai-(—l/2,u), a[_{—l/2,u) 
and a 3+ (l/2,u), a\ + (l/2,uu) do not, if the plate is surrounded by (absorbing) matter. 
When the plate is surrounded by vacuum, then the characteristic transformation and 
absorption matrices T(u) and A (a;) fulfill the matrix relation 11 

T(cj)T + (cj) + A(cj)A + (cj) =1 (5.13) 

(I, unit matrix), and from Eq. (5.8) it then follows that the output amplitude oper- 
ators also satisfy bosonic commutation relations. In fact, it can be shown that this 
is not only true at the very input and output ports of the plate but in the whole 
half-spaces on the left and the right. In this case both the input amplitude operators 
and the output amplitude operators can be regarded as being bosonic operators as- 
sociated with ordinary monochromatic incoming and outgoing modes, respectively, 
and the matrices T(cj) and A (a;) read as [n{uj) = 112(00)] 

T n (u) = T 22 {uo) = e-^ l/c r{uj) [l - t 1 (cu)e 2m ^/ c ^(cu)t 2 (cu)] , (5.14) 
T l2 {u) = T 2l {uo) = e- iul/ %{u)e m ^ l/c ${uj)t 2 {uj), (5.15) 

x [1 - e in ^ l/c r{u)} , (5.16) 



11 T 



i Note that a unitary transformation together with rescaling can be applied to the output ampli- 
tude operators such that the new operators are bosonic and thus Eq. (5.13) can be assumed to be 
valid, without loss of generality. 
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A 12 (u>) = -A 22 (u) = ^ ni {uj)n K {uj)e-^ 2c \{uj)d{uo) 



x 



y/X-(l,u) [l + e in( ^ l/c r(u)] . 



Here, 



and 



r{uj) 



1 — n(u) 
1 + n(uj) 



t 2 (u) = 



2n(u) 



(5.17) 



(5.18) 



(5.19) 



1 + n{uj) 1 ' 1 + n(u) 

are the interface reflection and transmission coefficients, respectively, and the factor 

#(u) = [l - r\u)e 2in( ^ l/c ] _1 (5.20) 

arises from multiple reflections inside the plate. 

Operator input-output relations of the type given in Eq. (5.8) are of course valid 
also for more complicated four-port devices such as multilayer plates. Obviously, the 
only difference consists in the actual expressions for the characteristic transformation 
and absorption matrices. For notational reasons it will be convenient to call the input 
amplitude operators di(u) [i.e., ai+(— 1/2, u)—>ai(u>), a 3 _(Z/2, a;)— > a 2 {uS)], the output 
amplitude operators bi(u) [i.e., ai-(—l/2,u>)—>bi(u), a 3+ (l/2, uj) — > b 2 {uj)], the device 
operators gi(u>) [i.e., g±(ou) — > gi{oj)] (see Fig. 2), and to introduce the definitions 




Figure 2: Scheme of a four-port device: Two incoming fields [photonic operators a\{uj) 
and a 2 (u)] are superimposed to produce two outgoing fields [photonic operators b\(uf) 
and b 2 {uj)], the §j(uj) being the operators of the relevant device excitations. 

a M= ( a f)Y b { u,)=( i f ) ), i( U )=f*^Y (5.21) 

The operator input-output relations (5.8) at a four-port device can then be written 
in the compact form of 

b(u) = T(cu)a(cj) + A(w)g(w). (5.22) 

In what follows we assume that the characteristic transformation and absorption 
matrices T(u;) and A (a;), respectively, obey the equation (5.13). 
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5.2 Quantum-state transformation 

The operator input-output relations (5.22) can be used to calculate various moments 
and correlations of the outgoing fields in a straightforward way [57, 59]. The output 
operators are expressed in terms of the input operators and the device operators 
which for themselves act on the quantum state the incoming field and the device are 
prepared in. Instead the photonic operators may be left unchanged but the quantum 
state is transformed. This equivalent procedure is suitable in view of problems, such as 
the determination of the entanglement of the outgoing fields, where knowledge of the 
quantum state of the outgoing field as a whole is required. Hence, we are interested in 
a unitary transformation that transforms the input-state density operator g in (i.e., the 
density operator of the quantum state the incoming field and the device are prepared 
in) into an output-state density operator g out according to 

Qout = Ug hl U\ (5.23) 

from which the density operator of the state the outgoing field is prepared in can be 
obtained by taking the trace with regard to the device variables. 

Let us shortly digress to lossless devices and assume, for a moment, that A(u) = 0. 
It is clear from Eq. (5.13) that in this case the characteristic transformation matrix 
T(u;) must be a unitary matrix and thus represents for each uj [and detT(c<j) = 
1] an element of the group SU(2) [60, 61, 62, 63, 64, 65]. For absorbing devices 
this can surely not be the case. Because of the coupling to the environment, we 
will definitely not be able to construct any unitary transformation that acts on the 
electromagnetic field operators alone. But we may look for one in the larger Hilbert 
space that comprises both the electromagnetic field and the device. Introducing some 
auxiliary field operators h(u) and defining the "four-vectors" 

we may extend the input-output relations (5.22) to the four-form of 

/3(w) = A(w)d(w), (5.25) 

where A(cj) is a unitary 4x4-matrix, hence A(u)A + (uj) —I. After separation of some 
phases from the matrices T(u) and A(u) and inclusion of them in the input operators, 
the matrix A(uj) can be regarded (for each a;) as an element of the group SU(4), and 
it can be expressed in terms of the matrices T(u) and A(u) [66] as 

{u) ~ \-s(u)c-\u)t(u) chs- 1 (cu)a(cu) ) ' (5,26) 



where 
and 



C(u) = y/T{u)T+(u) (5.27) 
S(u) = y/A{u)A+{u) (5.28) 
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are commuting positive Hermitian 2 x 2-matrices. Note that C 2 (uj) + S 2 (cj) = I. 

The matrix transformation in Eq. (5.25) may be realized also as a unitary operator 
transformation 

/9(w) = U ] a{u)U = A(u)a(u), (5.29) 

where 

U = exp!^-iJ^ <fa[a\u)] T &(u)a(u)^, (5.30) 

with the 4 x 4- matrix $(cu) being defined according to 

exp[-i*(w)] = A(w). (5.31) 

Obviously, [/ is just the unitary operator that transforms g in into g out in Eq. (5.23). 
Since the input density operator can be regarded as being an operator functional of 
6c(cu) and ct\ui), Qi n = Qi n [ac(ui), cJ(u)], from Eqs. (5.23) and (5.29) it follows that 
the transformed density operator can be given by 

g out = g- m [A+HaH, A T (w)a\uj\ . (5.32) 

Projecting g out onto the Hilbert space of the radiation field then yields the density 
operator of the outgoing fields 

eS = I & CD) {an[ A+ H & H» AT H« t ( w )]}» ( 5 - 33 ) 

where Tr*- -* means the trace with respect to the device. 12 

It is often useful and illustrative to describe quantum states in terms of phase- 
space functions such as the familiar s-parametrized phase-space functions (see, e.g., 
[68, 69]). From Eq. (5.32) it follows that the s-parametrized phase-space functional 
Pout[ot(w)', s] that corresponds to g ont is simply given by 13 

P out [a(u);s] = P iQ [A+(u)a(uj);s} , (5.34) 

so that the phase-space functional of the outgoing radiation reads as 

P ( 2[aH;s] = J VgP out [oc(u);s} = J Vg P in [A + (u)cx(u); s] , (5.35) 

where the functional integration (notation Vg) is taken over the continua of the com- 
plex phase-space variables gi(ui) and g2{oj) of the dielectric device. In Eq. (5.34) we 
have used the fact that application of the unitary transformation under consideration 

12 The similarity to the usual open-systems approach to dissipation is not accidental. However, 
master (or related) equations, to which an open-systems theory would lead [67], are not required 
here, because the action of the environment (i.e., the device) is explicitly known. For chosen quantum 
state the device is initially prepared in, all the necessary information is contained in the characteristic 
transformation and absorption matrices T(w) and A(w), which in turn are determined by the Green 
function of the phenomenological Maxwell equations of the classical problem. 

13 Since uj is continuous, P; n [a(o;);s] and P ou t[a.(uj); s] are functional rather than functions. 
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implies preservation of operator ordering, i.e., the annihilation and creation operators 
are not mixed by the quantum-state transformation. It should be pointed out that 
this is not the case for amplifying devices [70]. 

For practical purposes a description of the incoming and outgoing radiation in 
terms of discrete (quasi-monochromatic) modes is frequently preferred to be used. For 
this, we divide the frequency axis into sufficiently small intervals 14 of mid-frequencies 
uj m and widths Au m and define the discrete photonic input operators 



/ dua{uj), (5.36) 

J Au) m 



m J All)., 

and the discrete photonic output operators (3 m accordingly. Then we assign to each 
pair of operators ct m and (3 m the input-output relation (5.25) with the 4x4-matrix 
A m = A(cu m ), and, according to Eq. (5.30), the unitary operator U then reads as 

U = Y[U m , (5.37) 

m 

where 



U m = exp 



-«(ay T $ m d m (5.38) 



with <I> m = $(dj m ). In particular, the functional integral in Eq. (5.35) becomes an 
ordinary multiple integral. 



5.2.1 Examples 

Let us restrict our attention to (quasi-)monochromatic fields, so that it is sufficient 
to consider only a single frequency component u m . Suppose the incoming field and 
the device are prepared in coherent states 15 

|^„) = |7) = exp( 7 T a t -7 + «)|0), 7 = , (5.39) 

with Cj and dj (j = l, 2), respectively, being the coherent-state amplitudes of the input 
fields and the device. Application of Eq. (5.32) yields \ip out ) = \"f'), with 7' = A7 in 
place of 7 in Eq. (5.39), and it follows, on applying Eq. (5.33), that the outgoing 
fields are prepared in coherent states, i.e., 

£iHc')<c'|, c' = Tc + Ad. (5.40) 

Thus, the coherent-state amplitudes c[ and c' 2 of the outgoing fields are not only 
determined by the characteristic transformation matrix T but also by the absorption 
matrix A via the coherent-state amplitudes of the device. 

Next let us consider the case where the field in one of the two input channels is 
prepared in an n-photon Fock state and the field in the other input channel and the 

14 The matrices T(ui) and A(w) must not change substantially over an interval. 
15 For notational convenience we omit the subscript m. 
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device are left in vacuum, i.e., \ip in ) = |nOOO). The Wigner function of the input state 
reads 16 

W in (c*) = (-l) n 0) L n (A\ ai \ 2 ) exp ( - 2|a| 2 ). (5.41) 

Applying Eqs. (5.34) and (5.35) and integrating over the phase space of one outgoing 
field, the Wigner function of the field in the jth output channel is derived to be [66] 

WS (a,) = £ (f) \T^\ 2k (1 " \T 3l \r~ k W k (a,). (5.42) 

Obviously, Eq. (5.42) does not only hold for the Wigner function but for any s- 
parametrized phase-space function, and hence the corresponding density operator 
reads in the Fock basis as 

^ = E(fc)mii*(i-i^iiT"*i*><*i- ( 5 - 43 ) 

fc=0 ^ ' 

Note that the quantum state the outgoing field is prepared in contains only Fock states 
up to the photon number n of the input Fock state. That is a direct consequence of 
the action of the compact group SU(4) which leaves the total number of quanta 
unchanged. Furthermore, the Fock state with the highest photon number (i.e., the 
input state) appears in the output state with a weight |Tji| 2n , that is, the probability 
of finding the same number of photons after transmission through (reflection at) an 
absorbing device decreases as |Tji| 2n . In fact, for a device in the ground state, as it 
is the case here, there is some classical reasoning explaining this result. One could 
imagine that each of the n photons "feels" the effect of a transmission (reflection) 
coefficient smaller than unity, which for n photons amounts to the nth power of the 
transmission (reflection) coefficient. This reasoning fails when the device is prepared 
in anything else then the ground state. 

Finally, let us assume that the field in one of the two input channels is prepared 
in a Schrodinger-cat state and the field in the other input channel and the device are 
left in vacuum, i.e., 

|^) = -L(| 7 ) + |- 7 »|000), (5.44) 

where | 7 ) is a coherent state, and N = 2[1 + exp(— 2| 7 | 2 )] the proper normalization 
factor. Using Eq. (5.40) with c\ = ± 7 , C2 = d\ = d>2 = 0, it is not difficult to derive the 
density operator of the field in the jth output channel. The result is 

flli = ^{It^-iXt^-iI + I - 7^) (-7^1 

+ (| 7 T J . 1 )(-7^ 1 ] + |-7^-i)(7^i|)exp[-2| 7 | 2 (l - \T 3l f)] }. (5.45) 



16 W in (a) - W n (a 1 )Wo(a 2 )W (gi)W a (g 2 ), where W k (a) = (2/7r)(-l) fc i fc (4| a | 2 ) C - 2 l a l 2 [L k (z), La- 
guerre polynomial] is the Wigner function of a /c-quanta Fock state (see, e.g., [71]). 
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Whereas the two peaks decay as (T^l 2 , the quantum interference, in contrast, decays 
exponentially as [T^l 2 exp[— 2|7| 2 (1 — | T^- 1 1 2 ) ] - The larger the mean number of photons 
(n) = |7| 2 tanh |7| 2 becomes, the faster is the decay. Let us consider, e.g., the transmit- 
ted Schrodinger cat (j = 2). From Eq. (5.13) it follows that the squares of the absolute 
values of the transmission coefficient IT21 1 2 , the reflection coefficient IX22 1 2 , and the 
absorption coefficients |A 21 | 2 and \A 2 2\ 2 are related to each other as 1 — |T 2 i| 2 = |T 22 | 2 
+ \A 21 \ 2 +\A 22 \ 2 , so that exp[-2| 7 | 2 (l - \T 3l \ 2 )] = exp[-2| 7 | 2 |T 22 | 2 ) exp[-2| 7 | 2 (| A 21 \ 2 
+ I ^4.22 1 2 ) ] - The terms exp ( — 2 1 7 1 2 1 T22 1 2 ) and exp[— 2| 7 | 2 (|A 21 | 2 + |A 22 | 2 )] then respec- 
tively describe the decoherence associated with the losses owing to reflection and 
absorption. 



5.2.2 Entanglement degradation 

Quantum information processing such as quantum teleportation and quantum cryp- 
tography [72, 73, 74, 75, 76] is essentially based on entanglement, which can be 
regarded as being the nonclassical contribution to the overall correlation between 
two parts of a system. In particular, when two fields that are prepared in nonclas- 
sical states are superimposed by a four-port device, then the two outgoing fields 
are prepared in an entangled state in general. Entanglement of a bipartite quantum 
state g is commonly quantified by measures which fulfill some basic requirements 
as non- negativity (being zero only for separable states), invariance under local uni- 
tary transformations of the subsystems, and non-increase under arbitrary positive 
trace-preserving maps [77, 78] . Further, the reduced von Neumann entropy should be 
recovered for pure states which in addition gives a normalization condition for the 
measure. So far, the distance E(g) of g to the set S of all separable quantum states 
a 17 which is measured by means of the relative entropy has been proven to be one 
measure satisfying all the given conditions [77]. Thus, 

E(g) =minTr[£(ln£-m<7)] . (5.46) 

Let us consider the entanglement degradation that is observed when two fields 
that are prepared in a Bell basis state 18 

|^) = -L(|01)±|10)) (5.47) 

or 

|$ ± ) = i=(|00)±|ll)) (5.48) 

are transmitted through absorbing four-port devices prepared in the ground state. 19 
Note that the entanglement of each of the Bell basis states (5.47) and (5.48) is equal 

17 A quantum state a of a system that consists of two subsystems A and B is called separable 
if ~J2iPi^i A ^ ® °i B \ where <t| A ' and represent quantum states of the subsystems, and 

T, i Pi = 1 - 

18 For the entanglement degradation of a two-mode squeezed vacuum, see [79]. 

19 For an application to photon tunneling through absorbing dielectric barriers, see [80]. 
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to In 2 (sometimes also named 1 bit). Applying Eqs. (5.32) and (5.33) to the input 
states |*0in) = |^ /± ) ® |0) D and |^ in ) = j^) (g) |0) D (|0) D , ground state of the two 
devices), we easily find that the outgoing fields are prepared in the quantum states 

£;(I^ ± )) = |[(2-|T 1 | 2 -|T 2 | 2 )|00)(00|] 

+ | (T 2 |01) ± Ti| 10)) (T 2 *(01| ± T*(10|) , (5.49) 

£l(l$ =t )) = |[(i-|T| 2 )(i-|T 2 | 2 )|oo)(oo| 

+ |T| 2 (1 - |T 2 | 2 ) |10)(10| + |T 2 | 2 (1 - |T| 2 ) |01)(01|] 

+\ (|00) ± 7\T 2 |11)) (<00| ± (TiT 2 )*(ll|) . (5.50) 

Here and in the following the notation (T) n = (T) 22 = Ri and (T) 12 = (T) 21 = T\ 
for the elements of the characteristic transformation matrix T; of the Zth four-port 
device is used (/ = 1, 2) [cf. Eqs. (5.14) and (5.15)]. The entanglement contained in the 

(V) -L- (V) 4- 

quantum states of the outgoing fields, ^ U {(|\1/ )) and f^utd^ ))> can be estimated 
by using the convexity property of the relative entropy. 20 In Eqs. (5.49) and (5.50) 
the output state is written as a sum of separable states and a single pure state. Since 
separable states have zero entanglement by definition, the entanglement of the whole 
state is bounded from above by the reduced von Neumann entropy of the respective 
pure states. If we assume equal transmission coefficients of the two devices (T = T 2 
= T), the bounds are given according to [70] 

E[£J t (\^})} <|T| 2 ln2 (5.51) 

and 

E[Q^(\^))] <£[(! + m 4 ) In (1 + |T| 4 ) - |T| 4 ln |T| 4 ] . (5.52) 

Suppose the two fields are transmitted through (equal) optical fibres with perfect 
input coupling (R — 0), so that the transmission coefficient T may be given by 

rp _ ^inwl/c _ e in K ul/c e ~l/L ^ gg-j 

with / and L — cj {n\Lo) being respectively the propagation length through the fibres 
and the absorption length of the fibres. Substituting of this expression into the in- 
equalities (5.51) and (5.52) yields the dependence on I of the bounds of entanglement. 
The exact dependence on I of the entanglement degradation calculated by applying 
Eq. (5.46) is shown in Fig. 3. One observes that the states |$ ± ) decay faster than 
the states |^/ ± ). Since the device has been left in the ground state, we can again use 
some classical reasoning to explain this behaviour. When the two fields are initially 

20 E(J2iPiQi) < ^2iPiE(gi), where J2iPi — 1 [81]- The partition into separable states and a single 
pure state is not unique. Is has been proven, however, that for a pair of spin-i parties there exists 
a unique (sometimes called optimal) decomposition such that the weight of the separable state is 
maximized [82] . The reduced von Neumann entropy of the extracted pure state with corresponding 
minimal weight is then exactly the amount of entanglement, hence the inequality reduces to an 
equality. 
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Figure 3: Entanglement degradation of Bell basis states \^) (full curve) and 
(dashed curve) after transmission through absorbing fibres as a function of the fibre 
length. 

prepared in a state l^), then only a single photon (i.e., either the photon from the 
first field or the photon from the second one) is effectively subject to absorption. By 
contrast, the two photons get effected simultaneously, if the initial state is a state 



6 Spontaneous decay 

Spontaneous emission is a prime example of the action of ground-state fluctuations 
on physically measurable processes. Einstein [83] already pointed out that, in order to 
obtain the Planck radiation law, a process as spontaneous emission must necessarily 
be included in the theory of atomic decay. Later on, the radiation properties of an 
excited atom located in free space have been a subject of many studies. In particular, 
the rate of spontaneous emission of an excited (two-level) atom in free space is given 
by (the famous Einstein ^-coefficient) 

Fu _ u ld 2 (61) 

37rfe c 3 ' 

where d is the absolute value of the matrix element d= (1\(1a\u) of the atomic dipole 
operator cIa [Eq. (6.6)] between the upper state \u) and the lower state \l), and ua 
is the corresponding atomic transition frequency. When the atom is surrounded by 
(dielectric) matter, then the ground state felt by the radiating atom is changed and 
thus the rate formula (6.1) must be corrected 21 [17, 84, 85, 86, 87, 88, 89, 90, 91, 92, 
93, 94, 95, 96, 97, 98]. 

21 In the strong-coupling regime the decay becomes non-exponential and cannot be described by 
a rate (Sec. 6.4). 
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6.1 Time dependence of the atom— field system 
6.1.1 Hamiltonian 

Let us assume that the surrounding matter can be regarded as being a dielectric 
of given complex permittivity. Applying the minimal-coupling scheme (Sec. 4.1), we 
may decompose the Hamiltonian [Eq. (4.1)] of the coupled system consisting of an 
atom and the medium-assisted electromagnetic field as 

H = H A + H M + H int , (6.2) 

where 

^ = J2^- a + 1 2 /d 3 rp A (r)y3 A (r) (6.3) 
is the Hamiltonian of the atom, 

/poo 

d 3 r J duhu)&(r,u)f(r,u) (6.4) 

is the Hamiltonian of the electromagnetic field and the dielectric matter, and 

#int = - V ^ Pa A(r Q ) + / d 3 rp A (r)y3 M (r) (6.5) 
^ m a J 

a 

is the interaction energy. Here we have omitted the A 2 term. In the electric-dipole 
approximation, the first term on the right-hand side of Eq. (6.5) simplifies to 

"E— PaMr a ) = -\[d A ,H A ]A(r A ), (6.6) 

a 

where 

d A = ^2q a r a (6.7) 

a 

is the atomic dipole operator. Restricting our attention to a two-level system, the 
atomic Hamiltonian H A reduces to 

H A = hu u \u)(u\ + hu>i\l)(l\ = ^hu A & z + const., (6.8) 

where a z — \u)(u\ — \l)(l\. Thus we may further simplify Eq. (6.6) to obtain, on applying 
the rotating wave approximation and recalling Eqs. (3.46) and (3.47), 

- V ^p Q A(r a ) = -^EH + \r A )d + H.c. (6.9) 

(d real), with a being defined by a — \ l)(u\. 

For a neutral atom, the atomic polarization defined by Eq. (4.14) reduces in the 
dipole approximation to 

P A (r) = d A 5(r-r A ), (6.10) 
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and thus application of Eq. (4.24) yields (E^ = — Pj^/eo) 

J d 3 rp A (r)y3 M (r) = -d A E4(r A ), (6.11) 

which for a two-level atom in the rotating wave approximation reads 

J d 3 rp A (r)^ M (r) = -a t E{J +) (r A )d + H.c. (6.12) 

Eqs. (6.9) and (6.12) yield 

H int = -a^\r A )d + H.c. (6.13) 

6.1.2 Heisenberg picture 

From the Hamiltonian given in Eq. (6.2) together with Eqs. (6.4), (6.8), and (6.13), 
the Heisenberg equations of motion read [95] 

^ = |«7tE«(r A )d + H.c., (6.14) 
= iu A dt + ^E$(r A )da x , (6.15) 

f (r, u) = -iuf(r, u;) + < 4J G*(r, r A , u)d a. (6.16) 

cr y Tinea 

Substituting in the expression for the electric field strength [Eq. (3.38) together with 
Eq. (3.35)] for f (r, uj) the formal solution of Eq. (6.16), we derive, on using the relation 
(3.28), 

EW(r,t) = EitLe(r,t) 

+ / diu — ImG(r,r A ,cj)d / dr e^^^r). (6.17) 

TT^O Jo C 2 J t , 

Substitution of this expression into Eqs. (6.14) and (6.15) then yields a system of 
integro-differential equations for the atomic quantities. 

In the Markov approximation the integro-differential equations reduce to 
Langevin-type differential equations. It is assumed that [after performing the u- 
integration in Eq. (6.17)] the time integral effectively runs over a small correlation 
time interval r c . As long as we require that t — t'^>r c , we may extend the lower limit 
of the r-integral to minus infinity with little error. Further we require that r c be small 
on a time scale on which the atomic system is changed owing to the coupling to the 
(medium-assisted) electromagnetic field. In this case, in the r-integral in Eq. (6.17) 
the slowly varying atomic quantity ct{t)q iu1at can be taken at time t and put in front 
of the integral, thus 

EW(r,f) = E«jr,t) + *(*)— / dcu — Im G(r, r A , cu)d CK - cu) (6.18) 
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[C(x) = ttS(x) + iVx 1 ]. Substitution of this expression into Eqs. (6.14) and (6.15) 
yields 

(6.19) 
(6.20) 



b z = -r(l + a z ) + [| atEL + L(r A )d + H.c. 
*t = [i( WA - M - |r]at + ^EtL(r A )d^, 



where 



2u\didj 



Im Gij(r A , r A , w A ) 



is the rate of spontaneous decay of the upper state and 

^ = Ml <p [°°duj ^ ImG ^( rA > rA ' w ) 



JO C UJ-UJa 

22 



(6.21) 



(6.22) 



is the contribution to the Lamb shift. Note that from Eq. (3.44) it follows that 
Eq. (6.21) can be given in the equivalent form of 



27T f°° f 

T = didj / dcu (0|4(r A ,cu)4 T (r A ,^ A )|0), 
n Jo 



(6.23) 



which exactly corresponds to Fermi's Golden Rule (see, e.g., [100]). 



6.1.3 Schrodinger picture 

The above given equations of motion do not only apply to the spontaneous emission 
but are also suitable for the study of the evolution of a two-level atom driven by an 
external medium-assisted electromagnetic field. If the atom is initially prepared in 
the upper state and there is no driving field, then the use of the wave equation 



(6.24) 



may be more appropriate for the study of the motion of the coupled atom-field system 
[101]. According to the Hamiltonian in Eq. (6.2) together with Eqs. (6.4), (6.8), and 
(6.13), the state vector can be expanded as 



|*(0> = C^e-^u) <g> |0) 

- /d 3 r / duCu(r,u,t)e- i( - u+u,l *\l) <g> \i,r,u), 



(6.25) 



where |0) is the vacuum state of the fundamental fields /;(r, u), and \i,r,u) is the 
state, where one of them is excited in a single-quantum Fock state. It is not difficult 
to prove that the probability amplitudes C u and Cu satisfy the differential equations 

C u (t) = h= I duj %r I (IT ./dir..,;; 



xG Jl (r A ,r,cu)Q t (r,cu,t)e- i (— ^ 



(6.26) 



22 



For the the overall (vacuum) Lamb shift, see, e.g., [99]. 
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C H (r, u, t) = -A= ^ v^(^)G;(r A , r, u) C u (t)e^^. (6.27) 

We now substitute the result of formal integration of Eq. (6.27) [Cu(r,uj,0) =0] into 
Eq. (6.26). Making use of the relationship (3.28), we obtain the integro-differential 
equation 

C u (t)= [ dt'K(t-t')C u (t'), (6.28) 
Jo 

with the kernel function 

A A f°° 2 

K (t-l/) = -P^- / da;=VlmG' ij (r A ,r A ,a;)e-^ A ^*-*'). (6.29) 
hiieo Jo c 

Taking the time integral of both sides of Eq. (6.28), we easily derive, on changing the 
order of integrations on the right-hand side, 



C u (t)= [ dt'K(t-t')C u (t') + l 
Jo 



(6.30) 



[C u (0) =1], where 

K(t - = ^ f 4 lmG ^ r ^ [e-(— _ il . (6.31) 
hTT€ Jo c 2 i{u - U A ) L J 

Note that in the Markov approximation the kernel in Eq. (6.31) simply becomes 

K{t-t') = -\T + i5u, (6.32) 

where Y and 5u are respectively given by Eqs. (6.21) and (6.22). 

The equation (6.30) is a well-known Volterra integral equation of the second kind. 
It is worth noting that the integro-differential equation (6.28) and the equivalent 
integral equation (6.30) apply to the spontaneous decay of an atom in the presence 
of an arbitrary configuration of dispersing and absorbing dielectric bodies. All the 
matter parameters that are relevant for the atomic evolution are contained, via the 
Green tensor, in the kernel functions (6.29) and (6.31). In particular when absorption 
is disregarded and the permittivity is regarded as being a real frequency- independent 
quantity (which of course can change with space), then the formalism yields the 
results of standard mode decomposition (see, e.g. [102, 103, 104]). 



6.2 Atom near a dielectric surface 

As a first example, let us consider the spontaneous decay of an excited atom near an 
absorbing planar dielectric surface. To be more specific, the distance z to the surface 
of the atom is assumed to be small compared to the atomic transition wavelength. For 
real permittivity, this configuration has been studied extensively in connection with 
Casimir and van der Waals forces (see, e.g., [85, 88] and references cited therein). 



6 SPONTANEOUS DECAY 



35 



Further, it can be regarded as being the basic configuration in scanning near-field 
optical microscopy (see, e.g., [105]) and related applications. 

In order to calculate the decay rate T, Eq. (6.21), the imaginary part of the Green 
tensor of two (infinite) half-spaces with a common interface [106, 107] is required. 
When the atom's half-space is the vacuum, then the Green tensor in this half-space 
can be decomposed into the vacuum Green tensor and a reflection term, i.e., 

G M (r, r', u) = G^(r, r', u) + R k i(r, r', w), (6.33) 

where the vacuum Green tensor G^(r, r', u) is obtained from Eq. (A. 34) for q = uj/c. 
From Eq. (A. 36) it is seen that ImG^"(r,r',w) =0, and thus we obtain, on applying 
Eq. (A.38), 

ImG>,r )W )^4i. (6.34) 
Combining Eqs. (6.21), (6.33), and (6.34), we obtain 

r = T + ^Imfl M (r A ,r A)WA ), (6.35) 

where r is given by Eq. (6.1) 

In the coincidence limit, the reflection term Rm(t, r, u) can be given in the form 
[39, 108] 

Rxx(z, Z, to) = R yy (z, Z, LU) 



8nq 



poo ■ poo h, p 2ipz 

\ dkkpe 2i ^ z r p (k) + — dk——r s (k), (6.36) 
Jo 8n J (j 

A poo 2if3z 

R zz (z, z,uj) = — dk k 3 —— r p (k), (6.37) 
Jo P 

where r p {k) and r s (k) are the usual Fresnel reflection coefficients for p- and s-polarized 
waves (q—u/c, (3— y 1 q 2 — k 2 ; the origin of the co-ordinates is on the interface). When 
the distance z of the atom from the surface is small compared to the wavelength, qz 
1, then the integrals in Eqs. (6.36) and (6.37) can be evaluated asymptotically to 
obtain 

7? 1 n2 ~ l | 1 ("- 1 ) 2 | jg (n~l)(2n-l) Q( ) , , 

zz ~ I67rq 2 z 3 n 2 + 1 + Mzn{n+1) + 12tt n(n + 1) + { h 1 ] 

R xx = R = \R zz - HO^l + o(z). (6.39) 

yy 2 16nz n 2 + 1 3tt n + 1 V ; V ; 

[e — e(u), n = y/e(u)]. Inserting Eqs. (6.38) and (6.39) into Eq. (6.35) yields (w = w A ) 
[109] 

r = To - ( 1 + | )(— V , l2 + O(z-i). (6.40) 

8 V d?) \u A zJ |eK) + l| 2 1 ; 1 ; 

It is worth noting that the leading term ~ z~ 3 in Eq. (6.40) exactly agrees with the 
result of the involved microscopic approach in [91]. It typically occurs for absorb- 
ing media and reflects the possibility of nonradiative decay, if the excited atom is 
sufficiently near the medium. 
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6.3 Real-cavity model of spontaneous decay in a dielectric 
medium 

Let us consider the situation when an excited atom is placed inside an empty micro- 
sphere (embedded in an otherwise homogeneous dielectric) with a radius much smaller 
than the wavelength of the atomic transition. This might serve as a model for spon- 
taneous decay of a (guest) atom in an absorbing dielectric. 23 From simple arguments 
based on the change of the mode density it is suggested that the spontaneous emis- 
sion rate inside a nonabsorbing dielectric should be modified according to T = nT , 
where n is the (real) refractive index of the medium and T is given by Eq. (6.1) 
[110, 111, 112, 113]. Here it is assumed that the local field the atom interacts with is 
exactly the same as the electromagnetic field in the continuous medium. In reality, 
the atom is located in a small region of free space, and the local field is different from 
the field in the continuous medium which leads to the introduction of the so-called 
local-field correction factor £, thus giving T = n£T . Different models have been used 
to compute £. Here we want to concentrate on the (Glauber-Lewenstein) real-cavity 
model, which for real n leads to [10] 

3n 2 \ 2 

' (6.41) 



2n 2 + 1 



Experiments suggest that this model is a good candidate for describing the decay of 
substitutional guest atoms different from the constituents of the dielectric [114, 115, 
116]. 

Now let us turn to the question what will appear for absorbing media. According 
to Eq. (6.21), the rate of spontaneous decay is proportional to the imaginary part of 
the Green tensor in the coincidence limit (i.e., the two spatial arguments are to be 
taken at the position of the atom). In the real-cavity model the Green tensor that 
enters the rate formula is the Green tensor for the system disturbed by a small free- 
space inhomogeneity. Obviously, it can again be given in a form like in Eq. (6.33), 
where for a sphere (with the centre as origin of co-ordinates) the reflection term 
R(r,r',uj) reads [117] 

e,o n=l m=0 v v ' v ' 

Cf HM >m (r,^) ®M >m (r',^) 
+C n A V)N >m (r, ^) ® N >m (r', ^) 



x 



(6.42) 



where M» nm (r,fc) and N« nm (r, fc) are the (even and odd) vector Debye potentials 
defined by 

M >m (r, fc) = V x [^ nm (r, k) r] , (6.43) 

23 Such a model applies as long as the atom cannot "resolve" the atomic structure of the surround- 
ing matter. 
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N >m (r, k) = iv x V x [V> m (r, k) r] , (6.44) 

with 

(cos \ 
sm {m<t>) (6-45) 

[j n (/cr), spherical Bessel function of the first kind; P™(cos0), associated Legendre 
polynomial; for the rather lengthy expressions of the generalized reflection coefficients 
C^ N \ see [117]]. 

In particular in the coincidence limit r — > r' — > only the TM-wave vector Debye 
potentials N« 10 (r,fc) and N«n(r,A) contribute to R k i(r,r',uj) and we find that 

i4*(r,r,u;)| r=0 = ^-C? {u)5 kh (6.46) 
where the generalized reflection coefficient Cf (u;) reads 
C7fV) = 

[i + g(n + 1) — ig 2 n — g 3 n 2 / (n + 1)] e %e 
sin g — g{cos g + in sin g) + ig 2 n cos g — £ 3 (cos ^ — in sin g)n 2 / {n 2 — 1) 

(6.47) 

[n— \/ e{u) and g— Ruj/c, with i? being the cavity radius]. Hence, when the atom is 
located at the centre of the sphere, we may insert Eq. (6.46) into Eq. (6.35) to obtain 
(for uj = wa) the decay rate in the form of 24 [95] 

r = r [i + Re cf(w A )] ■ ( 6 - 48 ) 

As long as the surrounding medium can be treated as a continuum this result 
is exact and valid for all admissible cavity sizes, without restriction to transition 
frequencies far from medium resonances. If the cavity radius is small compared to the 
wavelength of the atomic transition (but large enough for the continuum approach), 
we may expand the generalized reflection coefficient Ci(ua), Eq. (6.47), in powers 
of g — Rua/ c to obtain 

, ,. , 9ei(w A ) 
i — i o 



+ 



|2e(cu A ) + 1| 2 \oo A R 



9eiK) [28 


s(u A )\ 2 + lQe R (io A ) + l] 


5\2e(u A ) + 1 


4 



LU A R 



9,,uU ' a) [4|5K)| 4 + 4e R (c A )|e(c A )| 2 + e^K)-e 2 K)] 



|2e(w A ) + l| 

H^ [4kM f+2£RMl | +0 (^). ( , 49) 



24 Note that when in the coincidence limit (at the position of the atom) the Green tensor has the 
form of Eq. (6.33), then the rate of the spontaneous decay has the form of Eq. (6.35). Since the 
vacuum Green tensor G^(r, r', uj) has no longitudinal imaginary part and the scattering component 
Rki(r,r' ,uj) is purely transverse, the longitudinal electromagnetic field does not contribute to the 
decay rate. 
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For ei(cu A ) = 0, i.e, when absorption is fully disregarded, in Eq. (6.49) only the third 
term survives and reproduces exactly the Glauber-Lewenstein local-field correction 
factor (6.41). For an absorbing medium terms proportional to R~ 3 and R~ l are ob- 
served, which can give rise to a strong dependence of the decay rate on the cavity 
radius. In particular, the leading term proportional to i?~ 3 can be regarded as cor- 
responding to nonradiative decay via dipole-dipole energy transfer from the atom to 
the medium. 

Examples of the dependence of rate of spontaneous decay on the atomic transition 
frequency are plotted in Fig. 4 for a model permittivity of Lorentz type 




Figure 4: The spontaneous decay rate T, Eq. (6.48), for the real-cavity model of an 
atom in an absorbing dielectric is shown as a function of the atomic transition fre- 
quency uja near a medium resonance for the model permittivity (6.50) (uj> = 0.46 cut, 
7 = 0.05 a>r [left figure] and 7 = 0.2c<j t [right figure]) and _R = 0.02Aa- 



s(uj) = 1 + - Uv 2 , . (6.50) 

It shows a band gap between the (transverse) resonance frequency ujt and the lon- 



gitudinal frequency = \J uj\ + u>p. One observes, that precisely in the gap region 
the spontaneous decay strongly increases due to non-radiative decay channels which 
become important in that frequency range. 



6.4 Cavity QED 

If the radius of the microcavity is not small compared with the wavelength of the 
atomic transition, the cavity can act as a resonator, and it is well known that the 
spontaneous decay of an excited atom can be strongly modified when it is placed in 
a microresonator (see, e.g., [118, 119]). There are typically two qualitatively different 
regimes: the weak coupling regime and the strong-coupling regime. In the weak cou- 
pling regime the Markov approximation applies and a monotonous exponential decay 
is observed, the decay rate being enhanced or reduced compared to the free-space 
value depending on whether the atomic transition frequency fits a cavity resonance 
or not. The strong-coupling regime, in contrast, is characterized by reversible Rabi 



6 SPONTANEOUS DECAY 



39 



oscillations where the energy of the initially excited atom is periodically exchanged 
between the atom and the radiation field. This usually requires that the emission is 
in resonance with a high-quality cavity mode. Recent progress in constructing certain 
types of microresonators such as microspheres has rendered it possible to approach 
the ultimate quality level determined by intrinsic material losses [120]. 

Let us consider the spontaneous decay of an excited atom placed inside a spherical 
three-layer structure (Fig. 5). The outer layer (r>R 1 ) and the inner layer (0 < r < R 2 ) 
are vacuum, whereas the middle layer (i?2 < r < R\), which plays the role of the 
resonator wall, is a dispersive and absorbing dielectric. In particular when a Lorentz- 
type dielectric is assumed whose permittivity is of the form (6.50), the wall would be 
perfectly reflecting in the band-gap zone, i.e., u-r <wa< w l, provided that absorption 
could be disregarded (7— >0). The Green tensor G(r, r\u) for a spherical three-layer 
geometry can be found in [117]. 




Figure 5: Scheme of the spherical microresonator. 



6.4.1 Weak coupling 



For r and r' inside the inner (vacuum) sphere, the Green tensor looks like the one 
for the two-layered spherical geometry considered in Sec. 6.3. It has again the form 
given in Eq. (6.33) together with Eqs. (6.42) - (6.45), but with different reflection 
coefficients (see [117]). Hence, when the atom is located at the centre of the 

inner sphere, then the rate of spontaneous decay can again be given in the form of 
Eq. (6.48). Obviously, for a sufficiently thick cavity wall, exp[— ini(ouA)(Ri — R2)w/c\ 
<Cl, the generalized reflection coefficient in Eq. (6.48) reduces to that one for the 
two-layered configuration [Eq. (6.47) with R—R2], and thus for a true microresonator, 
i?2<^A/c^> 1, the decay rate becomes [101] 



T ~ T Re 



= F 



n(u)A) — «tan(i?2^A/c) 
1 — iniuJx) tan(_R 2 ^A/ c ) 

n R (u A )[l + tan 2 (i?2a; A /c)] 



[1 + n\(ujA) tan(i? 2 ^A/c)] 2 + ^|(^a) tan 2 (i? 2 ^A/c) ' 



(6.51) 



From Fig. 6 it is seen that [for the model permittivity (6.50)] the rate of spontaneous 
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0.9 1 , 1.1 

co/co x 



Figure 6: The spontaneous decay rate T/Tq, Eq. (6.48), of an atom in the microres- 
onator in Fig. 5 is shown for R 2 = 30At, R\ — R2 — At, = 0.5o;t, and 7 = 10~ 2 wt- 
The curves in the inset correspond to 7 = 10~ 2 cu T (solid line), 7 = 2 x 10~ 2 u; T (dashed 
line), and 7 = 5 x 10 _2 c<j t (dotted line). 



decay sensitively depends on the transition frequency. Narrow-band enhancement 
of spontaneous decay (r/r > 1) alternates with broadband inhibition (r/T < 1). 
The frequencies where the maxima of enhancement are observed correspond to the 
resonance frequencies of the cavity. Within the band gap the heights and widths of the 
frequency intervals in which spontaneous decay is feasible are essentially determined 
by the material losses. Outside the band-gap zone the change of the decay rate is 
less pronounced, because of the relatively large input-output coupling, the (small) 
material losses being of secondary importance. 



6.4.2 Strong coupling 

The strength of the coupling between the atom and the electromagnetic field increases 
when the frequency of the atomic transition frequency approaches a cavity resonance 
frequency. In order to gain insight into the solution of the integral equation (6.30) for 
the strong-coupling regime, let us consider the limiting case of only a single cavity 
resonance (frequency uoq) being involved in the atom-field interaction. In this case 
we may approximate, on using Eqs. (6.21) and (6.48), the kernel function (6.31) by 



= -IrcSuce-^-^-^e- 5 ^-^, (6.52) 
and thus the integral equation (6.30) corresponds to the differential equation [101] 
C u (t) + [i(u c - cu A ) + 5u c ) C u (t) + \T c 5oo c C u {t) = 0. (6.53) 



where T c is the decay rate at the cavity resonance [i.e., ujx = ojc in the rate formulas 
(6.48) and (6.51], and 8loc is the width of the cavity resonance. For small absorption, 
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7 <C uj t , up, ujp/cuT, the resonance lines in the band-gap zone can be regarded as 
being Lorentzians, and 

8«>c = (6.54) 

Equation (6.53) reveals that (under the assumptions made) the upper-state proba- 
bility amplitude of the atom obeys the equation of motion for a damped harmonic 
oscillator. In the strong-coupling regime, where uja = ujc and Tc^Suc, damped Rabi 
oscillations are observed: 25 



\C u (t)\ : 



cos 2 (fit/2), 



(6.55) 



where the Rabi frequency reads 

ft = ^2Y c 5uj c . (6.56) 
Typical examples for the time evolution of the upper-state occupation probability 




Figure 7: The temporal evolution of the occupation probability |C u (t)| 2 of the upper 
atomic state is shown for R 2 = 30At, R\ — R 2 = At, = 0.5wt, = 1.046448o;t, 
r A T /(2c) = 10" 6 , and 7 = 10" 4 wt (solid line), 7 = 5 x 10^ 4 cj t (dashed line), 7 = 
1CT 3 u;t (dotted line). For comparison, the exponential decay in free space is shown 
(dashed-dotted line). 

are shown in Fig. 7. The curves are the exact (numerical) solutions of the integral 
equation (6.30) [together with the kernel function (6.31)] for the model permittivity 
(6.50) The figure shows that with increasing value of the intrinsic damping constant 
7 of the material of the cavity wall the Rabi oscillations become less pronounced. 
Physically, larger 7 means larger absorption probability of the emitted photon by the 
cavity wall and thus reduced probability of atom-field energy interchange. 

25 Note that in the opposite case where Tc^Suic the solution of Eq. (6.53) is \C u (t)\ 2 = c~ Tct for 

LO A =U C . 
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7 Extensions to other media 

So far we have focused our attention to quantum electrodynamics in absorbing, 
isotropic, non-magnetic, linear media. The quantization scheme outlined in Sec. 3 
can, of course, also be generalized to other media. The starting point is the basic 
formula (3.35) relating the operator of the electric field strength to the dynamical 
variables of the system composed of the electromagnetic field and the medium. Re- 
calling Eqs. (3.16) and (3.31), we may write Eq. (3.35) as 

E(r,u)=ic^o |dVG(r,r»j N (r», (7.1) 

where 



he 

l N ( r ^) = U V -^£i(r,w)f(r,w). (7.2) 
In what follows we briefly comment on possible extensions. 

7.1 Amplifying media 

Let us consider a dielectric which in some space and frequency regions is amplifying, 
but in other ones absorbing. Amplification can be described by a complex permittivity 
e(r,u) which exhibits the familiar properties as in the case of absorption, most im- 
portantly, it still fulfills the Kramers-Kronig relations (3.19) and (3.20), except that 
the imaginary part is negative, £i(r, u) < 0. It can be shown that the quantization 
scheme in Sec. 3 also applies to amplifying media if Eq. (7.2) is modified according 
to [40] 

j N (r, u) = u^\ £i (t,u)\ [e(ei)f (r, u) + 0(-£i)f t (r, co)] (7.3) 

[0(x), Heaviside step function], which reflects the well-known fact that amplification 
requires the roles of the noise annihilation and creation operators to be exchanged 
[121, 122, 123, 124, 125]. Substitution of this expression into Eq. (7.1) yields the 
representation of the (w-components of the) operator of the electric field strength in 
terms of the fundamental fields, which replaces Eq. (3.35). All the other formulas in 
Sec. 3 remain valid. In particular, the proof of the commutation relations (3.41) and 
(3.42) can be given in the same way as in Appendix B for absorbing media. 

7.2 Anisotropic media 

In anisotropic dielectrics the scalar e(r,u) has to be replaced by a tensor e(r, cu), 
which is symmetric for reciprocal media (see [58]), 

eij(r,u) = £ji(r,w). (7.4) 

Thus the Green tensor has to be determined from the partial differential equation 

V x V x G(r, r, lu) - — e(r, u)G(r, r', u) = S(r - r'), (7.5) 

cr 
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which replaces Eq. (3.24). The solution for a uniaxial bulk dielectric is given in 
Eq. (A. 39). The symmetry property (7.4) ensures that the reciprocity relation (3.27) 
is valid, and the generalization of Eq. (3.28) reads 

f uj 2 

j d 3 s — e m (s, uj)G ik {s, r,uj)G* jl (s,r',uj) = ImGy(r, r',u). (7.6) 
The extension of Eq. (7.2) to anisotropic, absorbing media then reads 26 



j N (r,c) = c^s I 1/2 (r,c)f(r,c), (7.7) 
where 

e\ /2 {r,uj) = 0{v,u)e^ l2 {v,u)0- l {v,u). (7.8) 

Here the orthogonal matrix 0(r, uj) transforms £i{r,u) into its (positive) diagonal 
form £j(r, uj), 27 

ei(r, uj) = _1 (r, u)ei(r, u)0(r, uj). (7.9) 

Substitution of Eq. (7.7) into Eq. (7.1) yields the (^-components of the) operator of 
the electric field strength expressed in terms of the fundamental fields, which replaces 
Eq. (3.35). The further procedure is same as outlined in Sec. 3, and it can again been 
shown that the commutation relations (3.41) and (3.42) are fulfilled. 

7.3 Magnetic media 

Let us now include a possible magnetic response of the matter in the quantization 
scheme, restricting our attention to isotropic media, so that the susceptibility /i(r, uj) 
is a scalar. We start from the constitutive relations, written in a non-standard form, 
which preserves the interpretation of the fields E and B as the fundamental electro- 
magnetic fields, i.e., 

D(r,w) = e E(r,uj) + P(r, uj), (7.10) 
H(r,w) = n B(r, uj) -M(r,w). (7.11) 
Here, the polarization is given according to Eq. (3.9), 

P(r, uj) = e [e(r, uj) - 1] E(r, uj) + P N (r, uj), (7.12) 

and the magnetization reads 

M(r, uj) = k [1 - k(t, uj)} B(r, uj) + M N (r, uj), (7.13) 

where k(t,u) = /i~ 1 (r, uj) and k q = /Iq 1 . It is not difficult to prove that the Green 
tensor now obeys the equation 



uj 2 

VxK(r,i<j)Vx — -e(r,uj) 

c 2 



G(t,s,u)=6(t-s). (7.14) 



26 Using the results of Sec. 7.1, it is not difficult to include in the scheme also amplification. 
27 Note that, depending on the type of anisotropy in the dielectric, the direction of the optical 
axes might be frequency dependent [126]. 
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Subsequently, the integral relation (3.28) has to be replaced by 

J d 3 s «i(s, u)d s n G ki (s, r, u) [d s k G* nj (s, r', u) - d s n G* kj ( S , r', uj)} 
f uj 2 

+ d 3 s — ei(s,u)G in (s,r,uj)G* n (s,r',uj) = lmGi j (r,r',uj), (7.15) 

which can be proved correct in analogy to the derivation given in Appendix A for 
non-magnetic matter. 

The contribution to j (r, uj) of the "magnetization" changes Eq. (7.2) to 28 



l N ( r ^) = ^£i(r,w)f e (r,w) + V x yj ^|«i(r, u)\ f m (r, uj), (7.16) 

with f e (r, uj) and f m (r, uj) being the bosonic vector fields associated with the electric 
and magnetic properties respectively. Here we have assumed that the medium is 
absorbing, so that /xi(r, uj)>0 [51] and thus «i(r, uj) = — /xi(r, uj)/\ji{r, uj)\ 2 < 0. 

Substituting for j N (r',a;) in Eq. (7.1) the result of Eq. (7.16), we obtain the (uo- 

components of the) operator of the electric field strength E(r, a;) expressed in terms 
of the dynamical variables of the system. The fields B(r, a;) and D(r, a;) can then be 
constructed in the same way as in Sec. 3, and by means of Eq. (7.11) the field H(r, uj) 
can be obtained accordingly. As before, the commutation relations (3.41) and (3.42) 
are fulfilled. 



7.4 Nonlinear media 

In general, the linear response of a dielectric medium to the electromagnetic field 
is only the first term in the series expansion of the polarization in powers of the 
electric field strength, and in strong fields nonlinear interactions play an important 
role. For example, the propagating of intense light in Kerr-type nonlinear media can 
lead to the formation of soliton-like pulses, the dispersion of the group velocity being 
compensated by the Kerr nonlinearity [127, 128, 129]. There have been a number of 
approaches to quantize the electromagnetic field in nonlinear dielectrics [24, 25, 26, 
27], but mostly absorption is not taken into account in a consistent way. 

A possible way to treat absorption is to include in the Hamiltonian (3.34) an 
appropriately chosen nonlinear term that may be thought of as being an integral over 
a nonlinear function of the medium polarization field P(r), and thus the extended 
Hamiltonian is given by 

H = y"d 3 r^ dujhujf j (r,Lo)i(r,Lo) + H NL (7.17) 

with 

Hnl = y"d 3 r/i NL [P(r)], (7.18) 
28 Note that Eq.(3.16) changes to j (r, lo) = -iujV_ H (r,uj) + V x _Al N (r,cj). 
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where 

P(r) = J duP(r, uj) + H.c. (7.19) 

with P(r, uj) being defined according to Eq. (3.9). 

The Heisenberg equation of motion of f(r, uj) then reads 

iM(r, uj) = [f (r, u)),H]= Tiu f (r, uj) + [f (r, uj), H nl ] , (7.20) 

which can be rewritten as cDf(r, uj) — ujf(r,uj), where 

u = idt + h-'H^, (7.21) 

with H^ L = [H NL ,0]. Hence the frequency-integrated (operator version of) 
Eq. (3.22) changes to [49] 

V x V x E(r) - K(r) E(r) = J duj iuj[i j N (r, uj) + H.c, (7.22) 

with K being the superoperator 

K(r) = 1 (i9 t + h^H^y e(r,idt + ^ 1 ^ L ). (7.23) 

A similar equation holds for the (transverse) vector potential, from which the B-field 
can be obtained. Obviously, the approach ensures that the commutation relations 
(3.41) and (3.42) are satisfied. 

For a linear medium, the superoperator K simply describes the effects of dispersion 
and absorption, but for a nonlinear medium it introduces via H^ L additional nonlinear 
coupling terms. An application of Eq. (7.22) to one-dimensional propagation of light 
in a Kerr- type nonlinear medium is considered in [49], and a simplified version has 
been used to study soliton propagation in absorbing optical fibres [130]. 
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A The Green tensor 
A.l Basic relations 

In Eq. (3.24), the Green tensor G(r, r', uo) as a function of r and r' can be regarded as 
being the matrix elements in the position basis of a (tensor-valued) Green operator 
G = G{uj) in an abstract 1-particle Hilbert space, 

G(r,r',u) = (r\G\r'). (A.l) 

The matrix elements of the position operator r are given by 

(r|r|r / ) = r<S(r-r / ), (A.2) 

and the matrix elements of the associated momentum operator p, 

[xi,pj\=i5ij, (A. 3) 

read 

(r|p|r') = - V5(r-r'). (A.4) 
Let H = H{uo) be the tensor-valued operator 

H = ip x ip x -q 2 I = p 2 I - p <g> p - q 2 I, (A. 5) 

where 



2 



g 2 = -5(r,c), (A.6) 



and I is the unit operator, 

(r\I\r') = 5(r-r'). (A.7) 
Equation (3.24) corresponds to the operator equation 

HG = I, (A.8) 

as can be easily seen. Writing down Eq. (A.8) in the position basis, 

J d 3 s (r| J H"|s)(s|G|r / ) = (r|/|r'), (A.9) 
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and using Eqs. (A.2) and (A.4), we are just left with Eq. (3.24). 

Now it is not difficult to prove Eq. (3.27). From Eq. (A. 8) it follows that the 
equation 

G = H 1 (A.fO) 
is valid, and thus we find, after multiplying it from the right by H, 

GH = J, (A.ll) 

which in the position basis reads 

J d 3 s (r|G|s)(s|if|r') = (r|/|r'). (A.12) 
Using again Eqs. (A.2) and (A.4), after some straightforward calculation we derive 



G ik (T,T f ,u)=5 ji 5(T f -r). (A.13) 



Comparing Eq. (A.13) with Eq. (3.25), we immediately see the Green tensor indeed 
obeys Eq. (3.27). 

In order to prove Eq. (3.28), we introduce operators O x defined by 

(0% = (0*)^Oj, (A.14) 
From Eq. (A. 8) it then follows that 

G t H t = I. (A.15) 

Multiplying Eq. (A. 8) from the left by G x and Eq. (A.15) from the right by G and 
subtracting the resulting equations from each other, we find 

G X {H -H t )G = G t -G, (A.16) 

which in Cartesian components reads 

& m (H mn - Hl m )G n] = G] t - G ir (A.17) 

From Eq. (A. 5) it is easily seen that 

H mn -Hl m = 6 mn (q 2X -q 2 ). (A. 18) 

Representing Eq. (A.17) [together with Eq. (A. 18)] in the position basis and recalling 
the relations 

(rlGjjO = (r'fe|r>* = G^(r', r, u) = G£(r, r', w ), (A.19) 
we eventually arrive at Eq. (3.28). 
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A. 1.1 Asymptotic behaviour 

The Green tensor Gy(r,r',w) is holomorphic in the upper half-plane of complex uj, 
because of the holomorphic behaviour of e(r,u). In order to study the behaviour of 
Gij(r, r', uj) for — > oo and — >(), we first introduce the tensor-valued projection 
operators 

i^i-P, iii = e®£ (a.20) 

p 

and decompose H, Eq. (A. 5), as 

H=(p 2 -q 2 )l ± -q 2 P. (A.21) 

Applying the Feshbach formula [131], we then may decompose the Green tensor 
operator G = H 1 , Eq. (A. 10), as 



G = H 1 = PCPhP^P 

+ [/ ± -/ ll (/ ll M ll )- 1 / ll M ± ]K[/ ± -/ ± M ll (/ ll M 11 )- 1 / 11 ], 



where 



(A.22) 
(A.23) 



k = [i^hi^ - ^hPCPhP^Phi^ \ 

It is not difficult to prove that 

pHp = -P q 2 P, Ph J x = -/"g 2 /^, (A.24) 

I L Hp = -I^P, I L HI L = I^p 2 - q 2 )l L . (A.25) 
Combining Eqs. (A.22) - (A.25) and using Eq. (A. 6), we obtain for G the expression 

^2 



G = ~—r PiPiP^P 



or 



+[/ ± - /"(/i^/ii)- 1 /"^]^^ - i^PCPiPy^P], 



where 



K 



oj 



UJ 



1^ I - '—a ] I x + =V /^/"(Jiif/ii)- 1 /"^^ 



(A.26) 
(A.27) 



[e = e(r,w)]. 

Now the desired limiting processes can be performed easily. For | c<^ | — > we find 
[iW= e W(T)=e(T,u = ti)\ 



lim ^-G = -PCPi (0) P)- 1 P 

|w|->0 C 2 



and 



hm g 2 G = -5( )/ll(/ll£W/")' 1 /". 



(A.28) 
(A.29) 
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For \lo\ -^oowe arrive at, on recalling that Pi 1 - = I ± P = and e(r, a;) — > 1 if 
oo, 



CO 



lim -^G = lim q G = —I. 



(A.30) 



\u)\^oo C |aj|^oo 

Obviously, the first term on the right hand of Eq. (A. 26) is the singular part of 
G for — >0. Performing in that term the Taylor expansion 

e(r,oo) = e (0) (r) + coe^(r) + ..., (A.31) 

where £^(f) is real and e^(r) is imaginary [see Eq. (3.21)], we find that 

ReGy(r,r',a;) ~ to' 2 (\co\ -> 0) (A.32) 

and 



ImGy(r, r', oo) ~ oo 1 — > 0). 



(A.33) 



A. 1.2 Isotropic bulk material 

For a bulk material of given permittivity e(oo) the Green tensor reads (p= r — r') 



where 



q(oo) = \fe{oo) — = [n R (oo) +ini(u 

c c 



Aiiq 2 (oo)p 



CO 



It can be split up into a longitudinal and a transverse part according to 



G ll (r,r' )W ) = 



1 



47rg 2 



(A.34) 
(A.35) 

(A.36) 



and 



G ± (r,T f ,u) 
1 



1 



3p<g>p\ 1 



+ 9 



+ 



Anq 2 
i 1 



P 2 / P 3 



3 \ p ® p 



K qp (qp) 2 (qpf) \qp (qp) 2 (qpf) p 2 
In particular, from Eq. (A. 37) it follows that 



LO 



ImG (r, r, oo) — lim Im G (r,r',u) — - — n R (oo)I 

r'^r one 



(A.37) 



(A.38) 
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A. 1.3 Uniaxial bulk material 



Following [132], the Green tensor for a homogeneous uniaxial dielectric can be given 
in the form of 



jqt(u>)p e 



£ e iqt(u)pe e iqt(uj)p 

e t A-np e 4np 

e iqt(uj)p e _g«<?t(w)p 



47riq t {uj) 



(p x c) <S> (p x c) 

(p x c) 2 
-c<g)c 2(p x c) <g> (p x c) 



(p x c) 2 



(p x cy 



(A.39) 



[p 2 = e c pe~ 1 p, g t = v /e^o;/c]. Here, the vector c is the unit vector parallel to the 
crystallographic axis of the medium, e c = e c (u) the (complex) permittivity in that 
direction, and e t = e t (uj) the (complex) permittivity transverse to it. Note that for 
e c = e t = e [and e^ 1 = e" 1 /] the Green tensor (A. 34) for isotropic bulk material is 
recovered. 



B Commutation relations 



Let us first consider the commutation relations of the electric field and the vector 
potential. Using Eqs. (3.38) and (3.46) together with Eqs. (3.35) and (3.47), recalling 
the commutation relations (3.32) and (3.33), and applying the integral relation (3.28), 
after some algebra we derive 



4(r),£,(r') 
A k (r),Mr') 



0. 



= o, 



e E k (r),Mr') 



— I dw^ImGtofr.s.w) 



7T 



(B.l) 
(B.2) 

^(s-r'). (B.3) 



In Eq. (B.3) the a;-integral can be performed by applying the rule 



oo poo 

duj ... = lim / da; ... . 



(B.4) 



Note that, according to Eq. (A. 33), lmGu(r,r' ,u) behaves like u^ 1 as uo approaches 
zero. Thus we may transform, on recalling Eq. (3.26), the cj-integral into a principal- 
part (V) integral, so that Eq.(B.3) reads 



e 4(r), ii(r') 



v rduu 2 

ft— I -G km {r,s,u) 



^(s-r'). 



(B.5) 



The evaluation of the ^-integral in Eq. (B.5) can be performed by means of 
contour-integral techniques. Since G fcm (r, s,cu) is a holomorphic function of u> in 
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the upper complex frequency half-plane with the asymptotic behaviour according 
to Eq. (A. 30), the u;-integrals can be calculated by contour integration along an in- 
finitely small half-circle |o;| —p, p^O, and an infinitely large half-circle \u\ = R, i?— > 
oo, in the upper complex half-plane 



j — < 



da; ... = lim 



da; . . . 



p^O / |w|= P 

aij>0 



lim / 

i?^oo J\u\=R 
u)j>0 



da; .... 



The definition of J- L (ID as given in Eq. (A. 20) implies that 

(r|i- L (ll)|r') = ^'1)^-0, 
and thus we may write, on applying Eq. (A.l), 

duj uj 2 



Ml 



duj uj 2 

jG km {r,s,uj) 

UJ c 1 



(r|G(o;)/ ± |r /v 



UJ c z 



Note that 



and 



$H(r) = -V <g> V 



1 



47r|r| 

S ± (r) = S(r)-dHr). 



From Eq. (A.28) it follows that 



lim %-GI 1 - = -/"(/'W'r 1 / 11 /- 







|w|->0 C 



(B.6) 

(B.7) 

(B.8) 

(B.9) 
(B.10) 

(B.11) 



(I"i ± = 0) and therefore the integral over the small half-circle vanishes. Finally, from 
Eq. (A. 30) we see that 

, ,2 

(B.12) 



or 



lim 



gi l = -i A 



Hence, 



j — t 



P — ^ (r|G(o;)^|r') 



a; cr 



2i [ 
Jo 



UJ 



duj - Im (r|G(a;)J ± |r / ) = iirS (y - r'), 



and the sought commutator reads 

£o4(r),A*(i-') 



(B.13) 



(B.14) 



Since the corresponding commutation relations for the displacement field [Eq. (3.40) 
together with Eqs. (3.37) and (3.35)] can derived in a quite similar way, we renounce 
the derivation here and immediately give the result 



Ac(r),A(r') 



0, 



(B.15) 
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Ac(r),i,(r') 



(B.16) 



This is the result we have expected, because the polarization P(r) = D(r) — £oE(r) 
is related to the degrees of freedom of the matter and it should therefore commute 
with the radiation field operators. 

Recalling the relations V x A = B, fl = -eoE" 1 , —V(p = E 11 , and P = D — e E 
and using the commutation relations (B.l), (B.2), (B.14), (B.15), and (B.16), we can 
then derive further commutation relations such as 



e E k (r), Bi(r') = -ih e Mm d^5(r - r), 



^^(r-r'), 

</Hr),4(r') 
0. 



0, 



(B.17) 

(B.18) 
(B.19) 
(B.20) 



It should be pointed out that the commutation relations given above are valid for 
both the medium-assisted electromagnetic field considered in Sec. 3 and the total 
electromagnetic field considered in Sec. 4. Needless to say that quantities of the 
medium-assisted electromagnetic field and quantities of the additional charged par- 
ticles commute. 



C Equations of motion 

In order to prove that the quantization scheme yields the correct Maxwell equations 
and the Newtonian equations of motion of the charged particles, let us consider, e.g. 
the Maxwell equation (4.8). Using the minimal-coupling Hamiltonian (4.1), we find 
that 

D(r) = — [D(r), H] = — dV / du hw [D(r), f f (r', u)f(r', u)] 

III III J Jq 

a 

+4 / dV [D(r),p A (r'^ A (rO]+ l|dV [D(r), p A (r^ M (r')] • (C.l) 

Recalling the definition (4.13) of the displacement field and applying the commutation 
relation (B.20), we may simplify Eq. (C.l) to 

£>(r) = iy dV J^fahu [D M (r),^,w)f(r>)] 

D M (r), [p Q - q a A(r a )] 



— V — 

a 

Eq \ ^ 1 

ih ^ 2m r 



7 <PA(r), [p a - g a A(f Q )]' 



(C.2) 
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Using Eqs. (3.39) and (3.40) together with Eqs. (3.36), (3.37), and (3.35) and recalling 
the basic commutation relations (3.32) and (3.33), we find that the first of the three 
commutators in Eq. (C.2) gives 

Iff 00 - 1 - 

— / dV / dcu Tiuj [D M (r), fV, u)f (r, u)] = — V x B(r) . (C.3) 
in- J Jo /-to 

In order to calculate the second commutator, we use the commutation relation (B.16) 
and obtain 

i E [ 6 M(r), [p« - g a A(r a )] 2 ] = -ft(r). (C.4) 

Finally, recalling the definitions of ^a(i") [Eq. (4.2)] and Pa(t) [Eq. (4.3)] and using 
the particle-operator commutation relations, we derive 



Eq \ 1 

-^E^7 V ^(r), [p Q -g a A(r a )]' 



-Ji(r). 



(C.5) 



In Eqs. (C.4) and (C.5), the transverse and longitudinal current densities j^M an d 
jl( r )) respectively, are defined by 



ji (ll) (r) = ^ ?Q f Q «5 ±(ll) (r-fa)+H.c., 



(C6) 



where the velocity operator r a of the ath particle is 

r a = \ [r a ,H] = — [p a - g a A(f a )] . 



(C7) 



Substituting Eqs. (C.3), (C.4), and (C.5) back into Eq. (C.2), we just arrive at the 
Maxwell equation (4.8). The Maxwell equation (4.6) and the Newtonian equations of 
motion (4.10) and (4.11) can be derived in a quite similar way. 



